Systematics of azimuthal asymmetries in 
heavy ion collisions in the lA GeV regime 



W. Reisdorf, ^' , Y. Leifels^, A. Andronic^, R. Averbeck^, 
V. Barret*^, Z. Basrak'^, N. Bastid^, M.L. Benabderrahmane ^, 
R. Caplar^, R Crochet ^ P. Dupieux^ M. Dzelalija'', 
Z. Fodor^, P. Gasik"^, Y. Grishkin^, O.N. Hartmann^ 
N. Herrmann K.D. Hildenbrand ^, B. Hong'^, T.I. Kang^''^, 
J. Kecskemeti^, Y.J. Kim^, M. Kirejczyk'^, M. Kis^'^, 
P. Koczon^, M. Korolija^, R. Kotte^ T. Kress ^, A. Lebedev^, 
X. Lopez T. Matulewicz™, M. Merschmeyer W. Neubert^, 

M. Petrovici^ K. Piaseckig'"\ F. Rami^ M.S. Ryu^, 
A. Schiittauf^ Z. Seres ^ B. Sikora'^, K.S. Sim\ V. Simion^ 

K. Siwek-Wilczynska™, V. Smolyankin ^, M. Stockmeier^, 
G. Stoicea^ Z. Tyminski'^, K. Wisniewski D. Wohlfarth^, 
Z.G. Xiao^'j, H.S. Xu^, I. Yushmanov\ A. Zhilin^ 

(FOPI Collaboration) 

^ GSI Helmholtzzentrum fiir Schwerionenforschung GmbH, Darmstadt, Germany 

^FOPI 

'^National Institute for Nuclear Physics and Engineering, Bucharest, Romania 

Central Research Institute for Physics, Budapest, Hungary 

^Glermont Universite, Universite Blaise Pascal, CNRS/IN2P3, Laboratoire de 
Physique Gorpusculaire, Clermont-Ferrand, France 

^ Institut fiir Strahlenphysik, Helmholtz-Zentrum Dresden-Rossendorf Dresden, 

Germany 

^Physikalisches Institut der Universitdt Heidelberg, Heidelberg, Germany 

^Institute for Theoretical and Experimental Physics, Moscow, Russia 

^Kurchatov Institute, Moscow, Russia 

^Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, China 

^Korea University, Seoul, South Korea 

Tnstitut Pluridisciplinaire Hubert Curien, IN2P3-CNRS, Universite de Strasbourg, 

Strasbourg, France 

^Institute of Experimental Physics, University of Warsaw, Poland 

^Rudjer Boskovic Institute, Zagreb, Croatia 



Preprint submitted to Elsevier 



15 December 2011 



Abstract 



Using the large acceptance apparatus FOPI, we study central and semi-central 
collisions in the reactions (energies in A GeV are given in parentheses): ^''Ca+^'^Ca 
(0.4, 0.6, 0.8, 1.0, 1.5, 1.93), ^^Ni+^^Ni (0.15, 0.25, 0.4), 96Ru+96ru (0.4, 1.0, 1.5), 
962r+962r (o.4, 1.0, 1.5), ^^Sxe+Csl (0.15, 0.25, 0.4), is^Au+^^^Au (0.09, 0.12, 
0.15, 0.25, 0.4, 0.6, 0.8, 1.0, 1.2, 1.5). The observables include directed and elliptic 
flow. The data are compared to earlier data where possible and to transport model 
simulations. A stiff nuclear equation of state is found to be incompatible with the 
data. Evidence for extra-repulsion of neutrons in compressed asymmetric matter is 
found. 
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1 Introduction 



Over the past nearly three decades the phenomenon of azimuthally asymmetric particle 
emission in high energy heavy ion reactions, called 'flow' in most of the relevant literature, 
has emerged as a general phenomenon observed over a very wide range of energies and 
systems. Ever since flow in heavy ion collisions was discovered [ip] it was believed that 
basic information on the nuclear equation of state, EoS, and other hot and compressed 
nuclear matter properties, such as viscosity, could be inferred from the data with use of 
hydrodynamical approaches [3|1|5] . This approach continues to be important in the highest 
energy regimes available today, see recent examples in PI7] . However, special measures 
have to be taken to handle both the initial and the final ('freeze-out') phases of the 
collision which by definition are far from the close-to-local equilibrium situation assumed 
in hydrodynamical theories. 

In the SIS energy range (O.lA — 2.0A GeV) which we cover here and which is expected 
to access nuclear matter up to 2.5 times saturation density, it was recognized early that 
microscopic transport theory that relaxes the local equilibrium assumption was necessary. 
Such dynamical theories suitable for heavy ion reactions were developed in the eighties [8] 
and continue to be improved piTOllTTllTMMllTSllTMTllTMQll^^ Since then 

many works confronting flow data with transport code simulations have been done. The 
situation towards the end of the nineties was reviewed in More recently a coura- 

geous attempt to present an EoS constrained by heavy ion flow data has found wide-spread 
attention [27]. However the authors state that they 'did not find a unique formulation of 
the EoS that reproduces all the data'. Other published experimental works on flow com- 
paring data with various transport codes come to similar conclusions not only 
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does one fail to reproduce all the data with a given code, there is unfortunately so far 
no really satisfactory agreement between the various theoretical approaches. Observables 
other than conventional flow parameterizations have been proposed with some success. 
In [3OII3T] the azimuthal modulations of radial flow were confronted with simulations us- 
ing different codes. While no definite conclusions were forwarded in [30] using a quantum 
molecular dynamics code [32], a clear preference for a 'soft' EoS was found ^31j using a 
Boltzmann-Uehling-Uhlenbeck (BUU) transport code [22||23] . A drawback of [31] was that 
conclusions were drawn from one code using one observable at one incident energy. A 
similar limitation holds for subthreshold kaon production [331IM] which shows sensitivity 
|35|l36] only in a narrow beam energy interval just below threshold and above sufficient 
detection probability. 

It is clear that a large and rather complete data base on all the observable aspects of 
heavy ion reactions spanning large energy ranges in sufficiently dense steps is needed in 
order to be able to perform a stringent testing of the quality of a given transport code. 
Only then can reliable and convincing conclusions on basic nuclear properties be derived. 
Recently our Collaboration, using the large acceptance apparatus FOPI [37|38] at the 
SIS accelerator in Darmstadt, Germany, has started a large effort to complement earlier, 
more specfic, works (see some references given in subsequent sections) by a systematic 
investigation encompassing 25 system-energies. 

A first paper written in this spirit [39] has offered a detailed overview of pion production. 
Pions are in this energy regime by far the most copiously created particle. This study 
contained also information on the pion's longitudinal and transverse rapidity ditributions, 
as well as the two forms of azimuthally asymmetric pionic flow: directed and elliptic 
(the deflnition will be recalled in the next section). In a second publication [ID] we have 
established a systematics for identifled light charged particles in the most central collisions, 
studying the freeze-out chemistry, the degree of stopping and radial flow. The present 
work, flnally, will present a systematics of directed and elliptic flow of these identifled 
light particles in the same, SIS, energy regime. 

In this energy regime the passing time and the expansion time are comparable. While this 
offers a possibility to obtain a handle on such basic properties as the sound propagation 
velocity, it also has the consequence that rather complex collision geometries result, as 
shown in early hydrodynamic simulations [S]. The presence of relatively cold 'spectator' 
matter at the time where the flreball is already expanding leads to cross talk between 
the various zones of very different local temperature. Also, quantum effects such as Fermi 
motion cannot be ignored as the initial rapidity gaps between projectile and target are not 
extremely large. The resulting flow patterns are rather different from those observed [H] 
for instance at the relativistic heavy ion collider RHIC, although, formally, one is using the 
same parameterizations in terms of Fourier expansion coefficients [12P3] , such as vi and f 2 
to be defined in section [2l To cope with the complexity of the flow in heavy ion collisions at 
SIS we will present an overview covering dependences on incident energy, system size and 
centrality. As some of the observables are strongly dependent on centrality, we make an 
effort to deflne centrality in a way that allows a good match to centralities in simulations. 
It will become clear furthermore that a separation of flow contributions by the mass and 
charge of the ejectiles covering differentially a large part of a sharply deflned phase space 
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section both in transverse and longitudinal dimensions is needed to be able to draw unique 
conclusions from comparison with transport model simulations. 

We also devote some effort to a more recent extension of this physics to include the isospin 
degree of freedom jl8|i44j by varying the system's isospin on one hand, and by measuring 
the flow differences of the isospin pair and ^He on the other hand (the tt^/tt^ pair was 
studied in [39]). 

In the sequel we start by briefly describing the FOPI apparatus used in this work, the 
centrality selection and by defining some of the most used terms in section [2l An overview 
of the experimental data on directed flow, section [3l and elliptic flow, section HI follows. 
An orientational confrontation of our data with transport model calculations will be given 
in section [5] and we will terminate with a summary of the main conclusions. 



2 Apparatus and data analysis 

The experiments were performed at the heavy ion accelerator SIS of GSI/Darmstadt using 
the large acceptance FOPI detector [37f38 j. The equipment and analysis methods in the 
present work largely correspond to those used in two of our earlier publications on pion 
emission |39] and central collision observables |10] in the same energy regime. We therefore 
confine ourselves to mentioning a few points that should be useful in making the present 
work reasonably self-contained. A total of 25 system-energies are analysed for this work 
(energies per nucleon, E/u, in GeV are given in parentheses): ^"^Ca+^^Ca (0.4, 0.,6 0.8, 
1.0, 1.5, 1.93), 58Ni+58Ni (0.15, 0.25), ^^Sxe+Csl (0.15, 0.25), ^^Ru+^^Ru (0.4, 1.0, 1.5), 
962r+96zr (0.4, 1.5), i^^Au+^^^Au (0.09, 0.12, 0.15, 0.25, 0.4, 0.6, 0.8, 1.0, 1.2, 1.5). 

Two setups were used for the low energy data {E/u < 0.4 GeV) and the high energy 
data. For the latter, particle tracking and energy loss determination were done using two 
drift chambers, the CDC (covering polar angles between 35° and 135°) and the Helitron 
(9° — 26°), both located inside a superconducting solenoid operated at a magnetic field of 
0.6 T. A set of scintillator arrays. Zero Degree Detector (1.2° -7°), Plastic Wall (7° -30°), 
and Barrel (42° — 120°), allowed us to measure the time of flight and also the energy loss. 

In the low energy run the Helitron drift chamber was replaced by a set of gas ionization 
chambers (Parabola) installed between the solenoid magnet (enclosing the CDC) and the 
Plastic Wall. More technical details can be found in Refs. [37|l5] . The combination of the 
two setups allowed us to identify pions and light charged particles (LCP, i.e. isotopes of 
hydrogen and helium) over a broad sector of phase space and to cover a large range of 
incident E/u since there is approximately a factor of twenty between the lowest and the 
highest energy. 

As in our earlier work, collision centrality selection was obtained by binning distributions 
of either the detected charge particle multiplicity, MUL, or the ratio of total transverse to 
longitudinal kinetic energies in the center-of-mass (cm.) system, ERAT. To avoid autocor- 
relations we have always excluded the particle of interest, for which we build up spectra, 
from the definition of ERAT. We estimate the impact parameter b from the measured 
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differential cross sections for tlie ERAT or the multiplicity distributions, using a geometri- 
cal sharp-cut approximation. More detailed discussions of the centrality selection methods 
used here can be found in Refs. [461I47] . We characterize the centrality by the interval of 
the scaled impact parameter 60 defined by 60 = b/bmax, taking bmax = 1.15{A](^ + A]^^) fm. 
This scaling is useful when comparing systems of different size. In this work we present 
data for four centrality bins: 60 < 0.15, 60 < 0.25, 0.25 < &o < 0.45 and 0.45 < 60 < 0.55. 
The total number of registered events was typically (0.3 — 1.5) x 10^ and was triggered by 
a multiplicity filter. Some minimum bias events were registered at a lower rate. 

Figures [1] and [2] give an idea of the yield and flow data covered by our apparatus for 
mass and charge identified LCP in the high energy run. In these plots the abscissa is the 
longitudinal (beam axis) rapidity y in the c.o.m. reference system and the ordinate is the 
transverse (spatial) component t of the four-velocity u, given by Ut = f3t'y. The 3-vector 
(3 is the velocity in units of the light velocity and 7 = 1/^1 — f3'^. Throughout we use 
scaled units yo = y/yp and Uto = Ut/up, with Up = f3p'yp, the index p referring to the 
incident projectile in the c.o.m.. In these units the initial target-projectile rapidity gap 
always extends from ?/o = — 1 to ?/o = 1- 

In Fig. [1] the upper left panel shows for lA GeV Au on Au collisions the proton yield 
distribution in the {yo vs Uto) plane. The centrality is 0.25 < 60 < 0.45. Two separate 
parts are visible: a predominantly forward {yo > 0) part obtained with use of the carefully 
matched Helitron and Plastic Wall and a mainly backward part using the matched CDC 
and Barrel detector systems. For particle types separated by charge only there is an addi- 
tional forward low-momentum sector available from the Zero Degree detector that is not 
plotted here, see Ref. [101 • While for the centrality selection the complete detector system 
is used, the flow analysis proceeds as follows (as illustrated in the various panels of Fig. [1]): 
To remove problem areas and edge effects visible in the upper left panel we apply rather 
restrictive sharp cuts as shown in the upper right panel. Taking advantage of the symme- 
try of the system we then apply reflection symmetry on the yo = axis, lower left panel 
and, finally, after some smoothing [18|^ we fill some of the small gaps that remained 
by two-dimensional interpolation. The resulting two-dimensional distribution covers most 
of the phase space except for low transverse momenta. It is in this phase space that we 
are able then to determine the flow fields describing the azimuthal emission dependences 
shown in Fig. [21 

Owing to collective flow phenomena, discovered experimentally in 1984 [T|2] , it is possible 
to reconstruct the reaction plane event-by-event and hence to study azimuthal correlations 
relative to that plane. We have used the transverse momentum method [19] including all 
particles identified outside the midrapidity interval \yo\ < 0.3 and excluding identified 
pions. Pions, in this context, were only found to be important at the highest SIS energy, 
but then contributed primarily to the reaction plane fluctuation. 

We use the well established parameterization [^2]H3] 

u = {-f,l3-f) ; Ut = I3a (1) 
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Fig. 1. Two-dimensional yield distributions (contours differ by a factor 1.5) of protons in collisions 
of Au+Au with lA GeV beam energy and centrality 0.25 < 5o < 0.45. The ordinate is the scaled 
momentum uto, the abscissa is the scaled rapidity i/q. Upper left panel: data after matching 
of various detectors and subtraction of pionic contributions (in the forward part). Upper right 
panel: data after applying sharp cuts to remove problem areas and edge effects. Lower left panel: 
data after symmetrisation of the forward/backward hemispheres. Lower right panel: data after 
bi-dimensional smoothening and closing of some gaps by interpolation. 

vo[l + 2vi cos(0) + 2f2 cos(20)] (2) 



Utdutdyd(j) 

vo = vo{y,ut); vi = vi{y,ut); V2 = V2{y,ut) 




(cos(20)) 



(3) 
(4) 



where is the azimuth with respect to the reaction plane and where angle brackets indicate 
averaging over events (of a specific class). The Fourier expansion is truncated, so that only 
three parameters, Vq, vi and V2, are used to describe the 'third dimension' for fixed intervals 
of rapidity and transverse momentum. The adequacy of this truncation for the lA GeV 
beam energy regime was already noted in Ref. [50] where this kind of Fourier analysis 
of azimuhal distributions was probably first used. Due to finite-number fluctuations the 
apparent reaction plane determined experimentally does not coincide event-wise with the 
true reaction plane, causing an underestimation of the deduced coefficients Vi and V2 which, 
however, can be corrected by dividing events into randomly chosen sub-events [49)151] : as 
explained in more detail in [15], we have used the method of Ollitrault [52] to achieve this. 
The correction factors (for the ERAT selection) are listed in a table that can be found in 
the Appendix. The finite resolution of the azimuth determination is also the prime reason 
why the measured higher Fourier components turn out to be rather small. 



Alternatively to the three Fourier coefficients, one can introduce the yields Qi, Q2, Q3 
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Q4 in the four azimuthal quadrants, of which only three are independent (on the average 
over many events) due to symmetry requirements 



Q2 = Qa (5) 

Qo = Qi + Q2 + Q3 + Qi (6) 

Q24=Q2 + g4 (7) 

More precisely, the flow axis is defined to be the median axis of quadrant Qi which extends 
from = —45° to 45° relative to the reaction plane located by definition at = 0. Thus, 
Qs is the 'antiflow' quadrant, while Q2,4 are the out-of-plane quadrants. 

The two equivalent triplets 

Vo, Vl, V2 < y Qo, Ql, (^24 

are related by 

Qo =2tt Vo ^ 6.28 Vq (8) 
Ql - Qs 2V2 



Qo TT 



Vl ^ 0.900 Vl (9) 



Q24 1 2 _ . . 

^ - - = — V2 - -0.637 V2 (10) 

Qo 2 TT 

These relations show that Vi is a dipole, while V2 is a quadrupole strength. Statistical (count 
rate) errors can be deduced with use of elementary algebra. The fact that Vi, as well as 
V2, are found to be non-zero, is generally called 'flow' in the literature and in particular 
the first Fourier coefficient is taken to be a measure of 'directed flow', while the second 
Fourier coefficient has been dubbed 'elliptic flow'. In the sequel we will occasionally also 
use instead of Vi{yo,Uto) (number weighted flow) the scaled transverse momentum (uto) 
weighted flow Uxoiuo, vto)- Most of the time we will plot elliptic flow with an inverted sign 
(i.e. — ^2) for reasons that will become clear later (see also the interpretation in terms of 
'quadrupole' strength given above). 

In the language of particle-particle correlations, which avoids the explicit use of a 'reaction 
plane' we are actually measuring the correlation of one chosen identified particle with all 
the other (measured) ejectiles. We are not trying here to separate so called 'non-flow' from 
'flow' contributions or to establish two-, four- or more particle-particle correlations. We 
take the point of view that 'non-flow' correlations ( such as resonance and cluster decays, 
which are not obviously disconnected from flow) should be properly taken into account in 
a realistic microscopic transport simulation. On the other hand multi-particle correlations 
are to a large degree revealed in our flow data for multi-nucleon clusters. 



In Fig. [2] we show the two-dimensional data deduced for the same reaction as in Fig. [T]for 
directed flow (left) and elliptic flow (right). The flat acceptance cut Uto > 0.4 preserves 
the backward/forward (a) symmetry of the apparent (vi) V2. In the low energy data we 
can use such a flat cut only if we require Uto > 0.8 for mass and charge identified particles 
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Vo Vo 

Fig. 2. Two-dimensional distributions of directed flow (left) and elliptic flow (right) for protons 
in collisions of Au+Au with lA GeV beam energy and centrality 0.25 < 60 < 0.45. The ordinate 
is the scaled momentum uto, the abscissa is the scaled rapidity i/q. The data are limited to 
uto > 0.4. In the left panel the colours delimit nine linear cuts varying from a maximum vi = 1.3 
to zero (at mid-rapidity) . For technical reasons the values in the backward hemisphere have been 
sign-inverted leading to a reflection symmetric (on the yo = axis) topology. In the right panel 
the colours for — i'2(y0) ^^to) again delimit nine cuts from a maximum around 0.31 reached for 
high momenta at mid-rapidity to a minimum -0.46 at |yo| = 1-5. The red colour marks values 
close to zero separating areas of differing sign. 



since the Helitron was not available yet for this part of our experimental campaigns. Small 
symmetry violations present in the original data have been corrected out. The quadrant 
method is well suited for this [39j • The corrections turn out to be critical for small isospin 
differences such as observed between and ^He. We shall therefore come back to this 
later in sections 13.51 and 14.61 



A few general remarks on accuracy will be made here. In the sequel, all our plotted 
data are given with errors that consist of a quadratic addition of systematic uncertain- 
ties and statistical fluctuations on the standard confidence level. If these errors are not 
visible they do not exceed the symbol sizes. Generallly the systematic errors of Vi and 
V2 are most important, but there are exceptions that can be seen when the errors tend 
to blow up for instance for high momenta, and for some of the Ca+Ca data which were 
obtained with lesser statistics. The estimates of systematic errors are based on observed 
forward-backward symmetry violations (a specialty of FOPI is that there is some overlap 
of measurements in the forward, resp. the backward hemisphers in the cm.), violations of 
the zero crossing of vi at mid-rapidity, discontinuities when switching sub-detectors, de- 
viations from the condition that the data for Z = 1 should be equal to the yield weighted 
sum of the separated isotopes (which require for identification additional sub-detectors) 
and incomplete isotopic separation (primarily /^He, see also the discussion in section 
3.5). These errors may vary somewhat with the system-energy, typically we have for vi 
errors of 0.01 (protons and deuterons), 0.015 (mass 3) and 0.02 (alpha's). For t>2 typical 
systematic errors are 0.005 (protons), 0.007 (deuterons) and 0.1 (mass > 2). 



8 



3 Directed flow 



Even under the contraints of symmetric heavy ion systems, the flow flelds Vi and V2 have 
complex multidimensional dependences: 



where E/u is the incident beam energy per mass unit, Agy^, Zgyg are the system mass and 
charge. A, Z is the ejectile composition. As a consequence a complete systematics encom- 
passes an enormous amount of information. It is out of question to present all this informa- 
tion in one readable paper: the chosen one-dimensional cuts through the flow topology are 
necessarily restrictive and perhaps even somewhat arbitrary. Nevertheless we believe that 
theoretical simulations that reproduce all of the presented results are likely to reproduce 
also the projections not shown here. Additional information can be found in some of our 
earlier publications on nucleonic flow at SIS energies [29|31f45f46f53f54f55|l56f57f58f59j . 

In the sequel we will flrst treat the centrality (60) dependence, then the ejectile mass 
[A) and the system size (Asys) dependence, the variation with incident energy {E/u) 
and, flnally, the isospin (Z/A), (Zsys/Asys) effects. Although direct comparisons with data 
obtained with different apparatus are exceedingly difficult owing to different apparatus 
biases and phase space cuts we will also present a few such comparisons. The next main 
section, on elliptic (^2) flow, will follow the same pattern. 

3.1 Centrality dependences 

In Fig.[3]we show a sample of our f i(?/o) data integrated over Uto, but with data for Uto < 0.4 
cut off (see section [2]). Only the forward hemisphere data are shown, the backward data 
are antisymmetric. It is known experimentally since some time [6OJ that directed flow is 
maximal at some intermediate centrality: for 60 = it must however converge to zero 
for symmetry reasons, but does not do so in experiment because of flnite Bq resolution. 
Nevertheless, below 60 ~ 0.5 the flgure shows that the variation with centrality is relatively 
modest as even the most central collisions (60 < 0.15 or 2.2% of the sharp radius cross 
section) still are associated with a signiflcant signal. The flnite resolution effect is only 
part of the explanation, since the other aspects of flow, V2, see section HI and stopping 
[IQ] are found to vary rather strongly with Bq. One may conjecture that the geometrical 
decrease of asymmetry with decreasing Bq is to some degree compensated by the increase 
of pressure expected from a more efficient stopping in more central collisions [40]. As we 
shall see in section \5\ sensitivity to the stiffness of the EoS is highest for the most central 
60 interval. 

Looking at the various panels of Fig. [3] a striking feature is the large similarity of the 
data in the shown energy interval (0.4 to 1.5^4 GeV) which covers the high energy part of 
the present study. This is in strong contrast to the behaviour at lower energy (see later). 



vi = vi{E/u,Asys,Z, 



bo,A,Z, yo,Uto) 
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Au+Au utO>0.4 




Fig. 3. Directed flow vi of protons in Au+Au collisions for different indicated centrality ranges 
and incident energies. The smooth curves are least squares fits to vi = wii •yo+'^ia "yo- Transverse 
4-velocities Uto below 0.4 are cut off. 
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The use of scaled rapidities is essential to see this. In general, for proton flow, which is 
strongly 'smoothened' by superimposed random motion, the data can be well described by 
two-parameter fits fi(yo) = ^ii ■ Vo + ^13 ■ Vo which are also plotted in the figure. These fits 
take care of smaller, more subtle, changes with incident energy that are visible on closer 
inspection of the f 1(2/0) data. 

Limiting ourselves to the information contained in the mid-rapidity slope parameter vu 
we obtain in a compact form data comparing the centrality (60) dependence of multi- 
plicity selected (MUL) with Ei^^T selected data. Fig. HI While the two types of selection 
widely seem equivalent for sufficiently large Bq, it is remarkable that for bo ^ the ERAT 
selections converge at all energies, except the lowest one (0.12^4 GeV, upper left panel) to 
a smaller value of vu. Rather than making a difference in bo resolution power responsi- 
ble for this, it is tempting to conclude that very high multiplicity events are less central 
than events with maximal stopping ERAT. The sorting according to MUL is based on the 
idea that the participant / spectator ratio is increased with increasing multiplicity because 
the hotter participants emit more particles. However, for very central collisions, where 
little spectator matter is present, this is no longer a sufficient impact parameter selection 
criterion. There stronger stopping (larger ERAT) is correlated with higher achieved den- 
sities. This in turn leads in the expansion phase to more achieved cooling, hence more 
clusterization. This picture is supported by simulations with transport codes varying the 
achieved density by varying the stiffness of the EoS, and by data varying with system size 
the achieved stopping jlO] . 

Our observation also shows that sorting data according to maximal multiplicity [HI] does 
not in the present energy regime lead to the same stopping as sorting them according to 
ERAT first, and then picking out a fixed cross section sample with the highest stopping 
|40] . The latter procedure is not an 'autocorrelation' as maintained in [61], but a legitimate 
procedure searching for maximal (not 'typical') stopping and hence, presumably, maximal 
achieved density. 

The low transverse momentum suppression that we make for charge and mass identified 
particles is resulting from our apparatus limitations (see section [2]). As this limitation 
is not present in our charge (but not mass) separated data, we can get an idea of the 
effect caused by transverse momentum cuts by looking at charge one and two data in our 
study. We show an example in Fig. |5] for Au on Au collisions at 0.8A GeV (centrality 
0.25 < 60 < 0.45) both with and without the Uto > 0.4 constraint. Here we choose the 
transverse momentum weighted flow u^o (the corresponding mid-rapidity slope is dubbed 
Uxoi) for reasons that will become clear when we try (subsection 13.61) to compare with 
older data in the literature. As can be seen, flow is signiflcantly larger when low transverse 
momenta are suppressed. The flow of particles that are only charge separated is difficult 
to interpret because it is a superposition from various isotopes with very different flow and 
hence is heavily influenced by the relative yields of these isotopes (see section |5l Fig. |36|) . 
Therefore in the next subsection we take a look at the flow of individual isotopes. 



^ Such slope parameters for the transverse momentum weighted flow were first introduced in 
Ref. [BD] 
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Fig. 4. Mid-rapidity slopes vu of directed flow of protons in Au+Au collisions for different 
indicated incident energies as function of average centrality. Centrality selection is done using 
either ERAT (closed red circles) or charged particle multiplicity MUL (open black circles). The 
smooth curves are least squares polynomial fits. Transverse 4-velocities uto below 0.4 are cut off, 
except for the two lowest energies (uppermost panels) where the cutoff is 0.8 (as indicated). 
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bo bg 

Fig. 5. Directed flow of Z = 1 (left) and Z = 1 (right) fragments in Au+Au collisions at 0.8A 
GeV. Shown is the centrality {ERAT) dependence of the mid-rapidity slope parameter n^oi with 
(full red circles) and without transverse 4-velocities itfo below 0.4 cut off. 

3.2 LCP dependences of directed flow 



We show the ejectile mass dependence of directed flow for centrality 0.25 < bo < 0.45 as 
a function of rapidity yo (integrated starting from Uto > 0.8 for the low energy runs up 
to E = 0.25yl GeV and Uto > 0.4 for the higher beam energies) in Figs. [6] and [7| and as 
function of transverse four- velocity Uto in a rapidity interval 0.4 < yo < 0.8 in Figs. [8] and 
[9l In these plots the mass three data are an average of the isospin pair ■^H/^He. The more 
subtle differences between and ^He will be discussed later. 

The data span the whole energy range measured for the Au on Au system (O.OQA GeV to 
1.5 A GeV) and hence contains a rather detailed information that represents a challenge 
for any theoretical attempts to reproduce all of them. One observes a gradual evolution 
of the shapes of the curves with the mass of the ejectiles: The variation with either yo or 
Uto is seen to be more pronounced as the mass is increased. The general increase of flow 
with mass is not new of course [28]I62|I63] and is expected even in pure hydrodynamics 
treatments [G^ where it is explained qualitatively in terms of random (thermal) motion, 
which is in velocity space largest for the lightest particles, on top of a common flow. 
This interpretation clearly does not then lead to (nucleon) number scaling which has been 
claimed to be proof of deconfinement of the constituents (constituent quarks at RHIC 
energies) [65]). It is obvious that the curves for the different masses (A=l to 4) are not 
equidistant along the ordinate and the different shapes as both the function of rapidity 
and transverse momentum (watch the location of the maxima) show that a perfect number 
scaling, even with transverse kinetic energy instead of Uto, would not emerge. We shall 
see later that the same conclusion also holds for elliptic flow. Of course due to cuts in 
the Uto {yo dependence) and in yo {uto dependence) a global scaling, if present, is not 
tested here. Qualitatively, we expect the purely hydrodynamic minded view to be an 
oversimplification of reality since the system is not governed by ideal hydrodynamics, see 
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Fig. 6. Directed flow vi of protons (black circles), deuterons (red squares), mass-three clus- 
ters (filled green triangles) and alpha-particles (pink asteriks) in Au+Au collisions at 0.09^4, 
0.12 A, 0.15 A, 0.25^4 GeV for centrality 0.25 < 6o < 0.45. The mass-three flow is obtained as an 
equal-weight average of and ^He. Low transverse 4-velocities uto < 0.8 are cut out. 

namely the incomplete stopping and the stopping hierarchy of clusters observed in [i0ll66] . 

The smooth curves plotted in Figs. |8]and|9]are just polynomial fits with maximally three 
parameters such as vi{uto) = vm ■ Uto + Vtu ■ u^q + vns ■ m^q- These fits can be extrapolated 
down to Uto = giving an idea of how the flow data might converge to zero at low transverse 
momenta. However, we find that the proton vi{uto) data at the lowest incident beam 
energies (Fig.|8]) seem to have an additional puzzling structure. One tentative explanation is 
that the bump is caused by protons emitted from heavier clusters, the proton then carrying 
a fraction of the flow of their parents. This would mean that late decays (evaporation) 
must be included in simulation codes. This could be of some non-academic importance for 
efforts to deduce isospin dependences of the nuclear EoS from neutron-proton difference 
flows [68|69] as it is the slow particles that are generally found to be sensitive to mean 
fields. Evaporation could also have a sizeable influence on the a data for the lower beam 
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Fig. 7. Directed flow vi of protons (black circles), deuterons (red squares), mass-three clusters 
(filled green triangles) and alpha-particles (pink asterisks) in Au-|-Au collisions at OAA, 0.6A, 
0.8A, l.OA, 1.2A, 1.5A GeV for centrality 0.25 < bo < 0.45. The mass-three flow is obtained as 
an equal-weight average of and ^He. Low transverse 4-velocities uto < 0.4 are cut out. 
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Fig. 8. Directed flow vi{uto) of protons (black circles), deuterons (red squares), mass three 
fragments (green full triangles) and alpha particles (pink asterisks) in Au+Au collisions at beam 
energies of 0.09^, 0.12^, 0.15A and 0.25^ GeV indicated in the various panels. The centrality 
is 0.25 < 6o < 0.45. The longitudinal rapidity is constrained to the interval 0.4 < i/q < 0.8. All 
smooth curves are fits to the data assuming vi{uto) = vtu • uto + Vti2 • u^q + vns • ufo- These fits 
are used to extrapolate to uto = 0. For the proton data only the high momentum parts have 
been taken into account for the fits. 



energies. 



3. 3 System size dependences of directed flow 



System size dependences of directed flow were flrst observed by the Plastic Ball group 
[67] . A sample of our present data is shown in Figs. [TO] to [131 The systems studied are 
indicated in the flgures and their captions. The incident energies are 0.4^ and l.bA GeV, 
the chosen centrality is 0.25 < 6o < 0.45. For identifled protons and deuterons both the 
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Fig. 9. Directed flow vi{uto) of protons (black circles), deuterons (red squares), mass three frag- 
ments (green full triangles) and alpha particles (pink asterisks) in Au+Au collisions at beam 
energies of 0.4^, 0.6A, 0.8A, l.OA, 1.2A and 1.5A GeV indicated in the various panels. The cen- 
trality is 0.25 < 6o < 0.45. The longitudinal rapidity is constrained to the interval 0.4 < yo < 0.8. 
All smooth curves are fits to the data assuming vi{uto) = vtu ■ uto + Vti2 • u'^q + I'tis • Uto- These 
fits are used to extrapolate to uto = 0. 



17 




Fig. 10. Left: Directed flow vi of protons (left) and deuterons (right) in collisions of Ca + 
"^^Ca (blue open triangles), ^^Ru + ^^Ru (black crosses) and ^^^Au+^^^Au (red dots) at incident 
energy of 0.4^4 GeV and centrality 0.25 < 6o < 0.45. The smooth curves are fits to vn yo + vi^ i/q. 
A cut uto > 0.4 is applied. 
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Fig. 11. Directed flow vi{uto) of protons (left) and deuterons (right) in collisions of ^*^Ca + ^^Ca 
(blue open squares), ^^Ru + ^^Ru (full black triangles) and ^^''Au+^^''Au (red dots) at incident 
energy of 0.4^4 GeV and centrality 0.25 < 6o < 0.45. The smooth curves are extrapolations from 
fits vtii ■ Uto + ^^412 • + vti3 ■ to the data in the range below uto = 1.2. A cut 0.4 < yo < 0.8 
is applied. In going from proton to deuteron flow there is a change of scale by a factor 1.44. 

rapidity dependences and the transverse momentum dependences are addressed. For the 
latter a rapidity interval 0.4 < i/q < 0.8 was selected. 



A few remarks can be made at this stage: 

1) Correlations between directed flow and stopping have been suggested earlier, [66]. By 
now extensive stopping data for the SIS energy regime exist [lO]. It will therefore be a 



18 




Fig. 12. Directed flow vi of protons (left) and deuterons (right) in collisions of Ca + Ca (blue 
open triangles), ^^Ru + ^^Ru (black crosses) and ^^^Au+^^''Au (red dots) at incident energy of 
1.5A GeV and centrality 0.25 < 60 < 0.45. The smooth curves are fits to vn yo + ^13 i/q + ^15 y^. 
A cut uto > 0.4 is applied. 
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Fig. 13. Directed flow vi{uto) of protons (left) and deuterons (right) in collisions of ^"^Ca + ^"^Ca 
(blue open triangles), ^^Ru + ^^Ru (full black triangles) and ^^^Au+^^^Au (red dots) at incident 
energy of 1.5A GeV and centrality 0.25 < 60 < 0.45. The smooth curves are extrapolations from 
fits vtu Uto + Vti2 u^Q + vti3 u^Q to the data in the range below uto = 1. A cut 0.4 < yo < 0.8 is 
applied. In going from proton to deuteron flow there is a change of scale by a factor 1.55. 

special challenge for future transport model simulations to try to reproduce simultaneously 
the system size dependences of stopping and of directed flow. 

2) The system size variations are of similar magnitude as predicted variations from mod- 
ifying the stiffness of the EoS (see section [5]): (partial) transparency can mockup a soft 
EoS. 

3) As before, we note the absence of simple-minded (nucleon) number scaling in the SIS 
energy regime for any system size. Looking at the data for many incident energies, we find 
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Fig. 14. Excitation function of directed flow vn of protons (left) and deuterons (right) in Au+Au 
collisions, vn was determined from a fit vi{yo) = vu ■ yo + V13 ■ in the range —1.1 < yo < 1.1. 
The various selections on the centrality 69 indicated. A cut Uto > 0.4 (0.8) is applied to the 
data with full (open) symbols. Note the ordinate scales differing by approximately a factor 1.5. 



ratios of 1.45 dz 0.10 for deuteron to proton flow. Since a switcfi to (transverse) kinetic 
energies for the abscissa would not modify the ordinates (fi), this conclusion holds also 
for the flow scaling as proposed by the RHIC community (which however was for elliptic 
flow, see later, section H]). 



3.4 Excitation functions for directed flow 



Excitation functions for the midrapidity slope of protons and deuterons in the Au on Au 
system are shown for two different centralities (60 < 0.25 and 0.25 < bo < 0.45) in figure 
[m In order to be able to joinup the higher energy data with the low energy data we apply 
a cut Uto > 0.8. A range limited to the high energies is also shown with a cut Uto > 0.4. 
The most striking feature is the steep rise up to OA A GeV, followed by a rather fiat 
behaviour beyond this energy. We stress that this shape is linked to the use of scale-free 
quantities, vi and yo- Repulsive pressure starts developing close to 0.1 A GeV, increasing 
first at a non-linear rate with incident energy. This regime is followed by a scenario where 
pressure augments at about the rate of the energy brought in causing little change in vu. 
The asymmetry caused by the sidepush is obviously larger if low transverse momenta are 
suppressed: this follows from the more detailed presentations shown before in subsections 
13.11 and 13. 2[ Comparing protons with deuterons, we note the change of the ordinate scale 
by approximately a factor 1.5, far from both one (no change) and two (nucleon number 
scaling) . 



20 



3.5 Directed flow and isospin 



The observed flow of particles is strongly dependent on tlieir mass as we fiave seen before. 
In our present study, confined to cliarged LCP's, tliere is only one pair, ^H/^He, that can 
potentially give information on isospin dependences of flow of ejectiles having the same 
mass. Another way to try to catch information on isospin is to vary the isospin of the 
systems. 

Before presenting some of our data relevant to the isospin degree of freedom, we feel it is 
necessary to point out in more detail some of our apparatus response affecting the results. 

Apparatus distortions influence the measured vi and V2 values and herefore have to be 
corrected for. The corrections turn out to be critical for small isospin differences such 
as observed between and ^He. For the symmetric systems that we study here, vi{yo) 
should be asymmetric with respect to midrapidity (?/o = 0), in particular fi(0) should be 
zero. As can be seen in Fig. [151 showing data obtained with use of the CDC (which covers 
predominantly the backward hemissphere) , vi{yo) does not cross the origin of the axes. 
We have checked that this mid-rapidity offset or 'f i-shift' was not due to deviations from 
azimuthal isotropy in the laboratory reference system which were below the 5% level and 
could be shown to produce, on the average over many events, only a negligible second order 
effect on the distributions with respect to the reaction plane. A systematic investigation 
of the wi-shift showed that it depended on particle type, centrality and system size in a way 
suggesting that it was correlated with the track density difference in the 'flow' quadrant 
Qi and the 'antiflow' quadrant Q^. While this could be simulated using our GEANT based 
implementation of the apparatus response, a sufficiently accurate quantitative reproduc- 
tion of the offset at mid-rapidity was not achieved. We therefore opted for an empirical 
method to correct for the distortion, using the very sensitive requirement of antisymmetry 
with respect to midrapidity. 

The quadrant method, defined in section [2l is well suited for this. We correct the Qi 
{i = 1, 2, 4, 3) by replacing them by Qi + CiQi. The coefficients q are fixed by three con- 
ditions. First, as we could reproduce total charges in central collisions to within 5% in 
25 system-energies |10] we require that after correction the total sum Qo is unchanged. 
Second, we assume a linear interpolation of the correction for (524, i-e. C24 = 0.5(ci + C3) 
and finally we demand f 1 = at ?/o = 0, i.e. Qi{l + ci)) = Qsi^l + C3) at i/q = 0. This is 
done separately for each particle, centrality and beam energy. A check of the Wi-shift as 
a function of pt (equivalently Uto) or of polar angle showed that satisfactory results were 
obtained assuming the distortions were proportional to pt in the range covered by the 
CDC, and independent of pt in the range covered by the Plastic Wall. 

In the upper right panel of Fig. [15] we show the effect of the correction for the mass three 
isospin pair. Clearly, the surviving isospin-effect is small. A slightly steeper slope of '^H, 
visible in the figure, was also found in the forward (Helitron-Plastic Wall) part of the data. 
Also, as we shall see in subsection 14.61 the f 2 data, which are affected by the corrections 
just described, also support more repulsion for than ^He. This gives us confidence in 
the consistency of the corrections. 
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Fig. 15. Comparison of the rapidity dependence of directed flow of and ^He in Au+Au 
collisions at l.OA GeV incident energy and centrality 0.25 < bo < 0,45 using the CDC filter. 
Upper left panel: uncorrected data using the CDC filter. Upper right panel: corrected data using 
the CDC filter. Lower left panel: uncorrected data using the full apparatus after symmetrization. 
Lower right panel: corrected data using the full apparatus after symmetrization. 



The two lower panels show the results when combining the data of all the subdetectors of 
FOPL In the left lower panels we ignore the f i-shift problem but symmetrize the data as 
outlined in section [2] with Fig. [H in the lower right panel the data from all the detectors 
have been first corrected to eliminate the Ui-shifts and subsequently were symmetrized. 
This is the procedure adopted for all the data shown in this work. 

The two ways of studying isospin signals are used in Fig. [16] comparing flow in the equal 
mass systems ^'^Ru + ^^Ru {N/Z = 1.182) and ^^Zr + ^^Zr {N/Z = 1.400) at an incident 
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Fig. 16. Centrality dependence of midrapidity slope vn of directed flow of the probes ^He 
(squares) and '^H (triangles) in the systems Ru+Ru and Zr+Zr at OAA GeV incident energy. In 
the upper two panels we compare the two probes in a fixed system as indicated in the panels, in 
the lower two panels we compare the same probe (^H left, ^He right) in two different systems. 

energy of OAA GeV and a centrality 0.25 < Bq < 0.45 with our usual flat (in yo) cut 
Uto > 0.4. If we look at the same ejectile, we find no convincing (within errors) difference 
for the mid-rapidity slopes Vn in the two systems (see the lower panels). Not shown here 
are data for all the other LCP and for a higher incident energy (1.5A GeV): these data lead 
to the same conclusion. If we compare the two ejectiles ^H/'^He in the same system (upper 
two panels), we see some evidence for a stronger flow of the more neutron rich cluster 
especially in the neutron richer system. In view of the systematic uncertainties definite 
and precise conclusions are premature, but the data are an encouragement for continued 
efforts to improve the statistics and the quality of such observables. 

Another way to search for isospin signals in directed flow is shown in Fig. [T7] where we 
plot for Au+Au {N/Z = 1.493) at six different incident energies (indicated in the various 
panels) the flow vi{uto) comparing (blue triangles) and '^He (red squares). Except for 
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Fig. 17. Directed flow vi{uto) of '^He (red open squares) and (blue open triangles) in col- 
lisions of Au+Au with centrality 0.25 < 60 < 0.45 and incident energies 0.12A, 0.25^4, OAA, 
l.OA, 1.2A and 1.5A GeV as indicated in the various panels. The longitudinal rapidity is 
constrained to the interval 0.4 < yo < 0.8. All smooth curves are fits to the data assuming 
viiuto) = ■Will ■ uto + Vti2 • Ufo + Vti3 ■ u^q. These fits are used to extrapolate to uto = 0. 
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the lowest shown energy, 0.12 A GeV (upper left panel), we always find the same pattern: 
a slightly higher fiow of in the range 0.4 < Uto < 0.8 with a tendency of the two data 
curves to merge at the higher values of Uto- The extrapolations below Uto = 0.4 have merely 
orientational value here. 

Closing this subsection, we can say that our data suggest that neutrons see a more repulsive 
mean field than protons in neutron rich compressed systems. This enhances somewhat the 
flow of clusters that contain more neutrons than protons. Of course it will be interesting 
to continue studies comparing single neutrons and protons with neutron sensitive devices 



3.6 Comparison to other data 



The comparison of flow data from different collaborations is a difficult task as both the 
centralities and the apparatus cuts or biases need to be reasonably well matched. Isotope 
separated LCP data have been published by the EOS Collaboration in 1995 [28j for the 
Au+Au system. The fiow quantity F defined in [28] corresponds to our midrapidity slope 
parameter u^oi ■ Pcm and hence has the dimension of a momentum (measured in GeV/c). 
We do not favour this 'half-scahng' (the scaled rapidity is used, but the momentum pro- 
jection is not scaled), but use it here to compare directly to the data in |28]. In Ref. |28] 
F was determined from a linear fit in the range —0.18 < yo < 0.29 and hence is techni- 
cally not strictly equal to our midrapidity slope which was determined in the larger range 
—1.1 < yo < 1.1 using a cubic (in yo) correction term in addition to the linear term, a 
procedure we found to be more robust against systematic errors. We note further that 
the recipee of centrality selection via reduced charged particle multiplicity bins used in 
[28] and earlier in [1] does not guarantee a fixed effective bo (or cross section) range as the 
incident energy is varied. This will influence the shape of the excitation function somewhat. 



For separated isotopes our measurements are affected by a transverse momentum con- 
straint, Uto > 0.4, which enhances flow as shown earlier using charge (but not mass) 
separated data that extend down to lower momenta. In Fig. [TS] we compare our flow data 
for the three hydrogen isotopes and for ^He with the EOS data [2H]- Clearly, due to the 
low-momentum cut our measurements lead to larger flow values, but come closer to the 
earlier data if we estimate the influence of the cut from the information given previously 
in Fig. [5l The remaining difference may be due to imperfect matching of centralities and 
a different data fltting range. Also, for lack of better knowledge, we have assumed that no 
matching apparatus fllter need to be applied before comparing. A more disturbing obser- 
vation is that our cut-corrected estimates for the isospin pair ^H/^He compare differently 
to the EOS data (see the lower two panels). For ^He, not shown here, the comparison is 
qualitatively similar to the '^He case. In contrast to what we tend to conclude from our 
data (see for example Fig. [T7j), the EOS data, presumably including low momenta, suggest 
a higher integrated F flow of ^He compared to ^H. It may be of some consequence that our 
reduction to zero uto cuts is only approximate since it was deduced from charge separated 
data only (Fig. [5]). 
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Fig. 18. Upper left: comparison of directed flow (midrapidity slopes) of protons in collisions 
of Au+Au (red, open circles) with centrality 0.25 < 60 < 0.45 (dubbed 'c2') with the data of 
Ref. [28] (black open squares). The latter are multiplicity selected in a bin dubbed 'M3+M4' 
(not further characterized in terms of cross sections or effective impact parameter). Upper right: 
same for deuterons. Lower left: same for tritons. Lower right: same for ^He. A cut uto > 0.4 was 
applied to the FOPI data (in red). The influence of this cut is estimated and the corresponding 
'rescaled' curves are plotted as well ( blue dots). The smooth curves are polynomial fits to guide 
the eye. 



4 Elliptic flow 



As announced earlier, we proceed to elliptic flow using a similar presentation of the various 
projections as was done for directed flow. But first we discuss briefly the question of the 
proper reference frame. 
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4-1 Question: to rotate or not to rotate 



One can rotate away the Vi component by an appropriate choice of the so-called 'flow' angle 
as new axis of symmetry instead of using the beam axis. This has been done in the early 
discovery days [50] and later also in other works [301I31II45] . The rotation generally leads 
to larger V2 components. While this is an instructive exercise, it does not add new data 
of course and is associated with a number of drawbacks. Technically, the Ollitrault [52] 
correction applied to V2 data is not designed for a rotated system. The most important 
drawback however, is the loss of rotational symmetry that the experimental equipment 
only has in the beam axis system. Due to imperfect An coverage this leads to incomplete 
(f) distributions in the rotated system and to a loss of othogonality of the various Fourier 
components. A more subtle detail is that the rotation ('flow'-) angle is not really a fixed 
global constant but is, besides being centrality dependent, strictly speaking also a particle, 
and even {yo,Uto) dependent quantity reflecting the multidimensional complexity of Vi. If 
one compares data with simulations where one wishes to convince oneself that both Vi and 
V2 are reproduced, it is an unnecessary complication to have to operate with rotations also 
in the simulation. For all these reasons we shall therefore stick to the beam axis system 
paying the price of a perhaps more complex interpretation in this system of reference. 
Once all the unrotated data are reproduced by the simulation, one can of course operate 
with rotations on the simulated events which have ideal in coverage, to try to extract 
possibly 'simple' interpretations. 



4-2 Centrality dependence of elliptic flow 

We show in Figs. [19] and [20] the rapidity dependence of proton elliptic flow for six different 
indicated incident energies spanning one order of magnitude: 0.15A to 1.5A GeV. In each 
panel the data for three centralities, bo < 0.25, 0.25 < 6o < 0.45 and 0.45 < 6o < 0.55 
are given. The shown energy regime is characterized by preferred out-of-plane emission 
at midrapidity (?/o = 0) 0: we plot (— ^2) here and in the sequel to associate this effect 
with a positive sign. We remind of the simple interpretation Q24/Q0 ~ 1/2 in terms of 
the quadrants. In the pioneering experimental papers [50]f70] the term 'squeeze-out' was 
coined describing the scenario in more visual language. Even earlier 'squeeze-out' was 
used in a theoretical paper [7T] with predictive power. The strong centrality dependence 
seen in the data suggests, alternatively, that we are observing an increasing shadowing 
due to the presence of spectator matter in the reaction plane in less central collisions. 
While this view point stresses the interaction between participants and spectators, the 
other viewpoint, squeeze-out, sees more the participant-participant interaction: already 
compressed ('shocked') matter being pushed out of the way by still incoming matter. The 
common feature of 'shadowing' and 'squeeze-out' is the fact that the expansion time scale 
and the passing time scale are comparable, a feature that is unique to the present energy 
regime, see also subsection 14.51 This creates sensitivity to properties tightly connected 
with the nuclear EoS, namely the sound propagation velocity (see section [5]). 



^ recall however the sign change at still lower energies close to the 'transition energy' [JS] 
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Fig. 19. Elliptic flow —V2 of protons in Au+Au collisions for different indicated energies and 
centrality ranges. The smooth curves are least squares fits to V2 = W20+f^22'yo+^24'yo- Transverse 
4-velocities uto below 0.8 are cut off. 

There are two other striking observations: 

a) The spectacular variation in the low energy regime, Fig. [191 followed by the relatively 
steady behaviour at the higher energy, Fig. [201 Note that in order to be able to follow the 
evolution in both regimes under similar constraints, we have used a common cut Uto > 0.8 
(upper panels in Fig. [2D1 the lower panels show data for a cut Uto > 0.4). 



b) The quasi-parabolic shapes leading to a sign change at large rapidities. The smooth 
curves are three parameter fits with even polynomials as indicated in the captions. Clearly, 
this observed shape is associated with using the beam axis as symmetry axis: we shall see, 
however, in section \5\ that the full rapidity range for ^2(2/0); not just the mid-rapidity 
part, is sensitive to the EoS allowing to make a more convincing and robust comparison to 
simulations than using just the value at yo = 0. In works such as [30], where a rotation of 
the symmetry axis was performed minimizing the cross correlations to directed flow, one 
finds 'squeeze-out' over a broader rapidity range, but at the cost of a rapidity dependent 
rotation ('flow') angle. As we shall see, we will be able to reproduce both ^1(^0) and V2{yo) 
with the same EoS. 

Additional information on the elliptic flow of protons is given in terms of its momentum 
dependence in the various panels of Figs. [211 and [221 for the same three centralities as 
before. A rapidity window |?/o| < 0.4 is taken, the smooth curves are three-parameter flts 
to V2{uto) = V22 ■ u^Q + V23 ■ u^o + f24 ■ uto i^owei energies) and to V2{uto) = V22 ■ u^q + V2i ■ 
ufo + ^26 ■ (higher energies). Again, the large differences between the low and the high 
energy regimes are evident. At the lowest three incident energies there is a momentum 
and centrality dependent sign change. At O.OOA GeV, and more weakly at 0.12y4 GeV, 
the most central collisions seem to display rather erratic fluctuations which, however, 
could beneflt from smaller error bars. Such 'critical' fluctuations could be resulting from 
a scenario where repulsive and attractive forces of almost equal strength on the average 
strongly compete. Further data from our earlier study of the complex behaviour in the 
vicinity of the transition energy are available ^5] . 
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Fig. 20. Elliptic flow —V2 of protons in Au+Au collisions for different indicated energies and 
centrality ranges. The smooth curves are least squares fits to V2 = f20 + V22 • Uq + ^^24 ■ Uq- In the 
upper three panels transverse 4-velocities Uto below 0.8 are cut off. The cut off is located at 0.4 
for the lower three panels. 

At the energies beyond 0.25^4 GeV a more regular pattern has emerged. The variation with 
centrahty is strong, stressing the need to match centralities very carefully when comparing 
different experiments or comparing to a simulation. One notices a gradual tendency of the 
shape changes to be more 'compact' as the energy is raised. One obvious reason for this is 
that, due to higher multiplicity, the phase space extension, in the scaled units used here, 
has to shrink since the total energy must be conserved. 



4-3 L CP-dependence of elliptic flow 



Moving to the cluster mass dependence of elliptic flow, we show in Fig. [23] a sample of 
the rapidity dependences for protons (black circles), deuterons (red squares) and alpha 
particles (pink asterisks) together with even polynomial fits (smooth curves). The chosen 
centrality is 0.25 < 60 < 0.45. We present the higher energy data (0.4 to l.bA GeV) 
with our usual constraint UtQ > 0.4. Again quasi-parabolic shapes are seen, but they are 
more compact in the rapidity direction as the cluster mass is increased. Whereas there is 
a clear mass hierarchy near midrapidity at the lower energy end (0.4^4 GeV, upper left 
panel), this hierachy is almost gone at 1.5^4 GeV (lower right panel), but the increasing 
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Fig. 21. Dependence of elliptic flow of protons on scaled moment uto for three indicated centrality 
bins in Au on Au collisions. The incident energies are 0.09^4 GeV (upper left), 0.12 A GeV (upper 
right), 0.15A GeV (lower left) and 0.25^ GeV (lower right). The rapidity interval is constrained 
to |yo| < 0.4. The smooth curves are 3-parameter fits to f2('Uto) = Vt22 • ufo + ^t23 • Uto + '"t24 ' ^to- 

compactness is preserved. This shows that limiting the flow information to mid-rapidity 
could be misleading when looking for possible interpretations. 

Elliptic flow data for the lower energies (O.OOA to OAA GeV) have been shown in an 
earlier publication of our Collaboration |15], but it is useful to complement the data given 
there with more details on the transverse momentum dependences, limited here by the 
longitudinal rapidity cut |?/o| < 0.4 and taken in the beam axis reference system. This is 
done in Figs. [24] and [251 As before for protons varying the centrality, we also find, at a 
fixed centrality, varying the particle mass, a fast evolution of the flow with incident energy. 
As observed in [IS], there is a sign change at some transverse four- velocity that depends 
on the particle type. The even polynomial fits (constrained to zero for Uto = 0) shown 
in the figures offer educated guesses of the flow behaviour down to zero momentum. At 
0.25^4 GeV it seems that a fully repulsive regime is reached. 
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Fig. 22. Dependence of elliptic flow of protons on scaled moment uto for three indicated centrality 
bins. Left: Au+Au at OAA, 0.8A, 1.2A GeV, right: at 0.6^, 1.0, 1.5 GeV. The smooth curves 
are 3-parameter fits to f 2(^*0) = 1^*22 • Uto + Vt2i ■ uto + '^t26 ■ u%. 
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Fig. 23. Elliptic flow —f 2(2/0) of protons (black circles), deuterons (red squares) and alpha par- 
ticles (pink asterisks) in Au+Au collisions for different indicated incident energies and centrality 
0.25 < 60 < 0.45. Transverse 4-velocities uto below 0.4 are cut off. The smooth curves are even 
polynomial fits. 
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Fig. 24. Left: elliptic flow —^2(^*0) of hydrogen isotopes for Au+Au at various indicated incident 
energies with centrality 0.25 < 60 < 0.45 and in the rapidity window |yo| < 0.4. The smooth 
curves are three-parameter polynomial fits including even (2,4,6) exponent terms. 



In Fig. [25] we illustrate again the sign change of f 2 at some (scaled) transverse 4- velocity 
Uto characteristic for the mass of the emitted cluster (upper left panel). It is possible 
to superimpose the flow data for the three particles if we introduce a rescaling of the 
abscissa by a constant factor, 0.56 for the proton data, 0.79 for the deuteron data, with 
the triton data left unchanged. (The ordinates (— ^2) are left unchanged). In the lower 
two panels of the same figure we show the effect of the same rescaling of the abscissa for 
the reconstructed contrained transverse rapidity, i/xmo, distributions |10] of the hydrogen 
isotopes. 'Constrained', denoted by the index 'm' for mid-rapidity, means that a cut |?/o| = 
lUzol < 0.5 has been applied. (Here also the ordinates have been renormalized to give 
the same integrated yields). This suggests that a random, perhaps thermal, motion is 
superimposed on a common flow profile at this incident energy. 



The transverse four- velocity dependences of elliptic flow, f2(wto), for the hydrogen isotopes 
at the higher incident energies (0.4^4 to 1.5^4 GeV) are shown in Figs. [26] and [271 again, for 
centrality 0.25 < 60 < 0.45. In Fig. [26]a midrapidity constraint \yo\ < 0.4 was applied, while 
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Fig. 25. Upper left: elliptic flow —V2{uto) of hydrogen isotopes for Au+Au at 0.15 A GeV with 
centrality 0.45 < 6o < 0.55 and in the rapidity window |?/o| < 0.4. The smooth curves are 
three-parameter polynomial fits including a quadratic, cubic and quartic term. Flow is seen 
to change sign for uto = 1.26 (protons p), 0.92 (deuterons d) and 0.73 (tritons, t triangles). 
Lower left: Transverse rapidity distributions [4ilj of the same particles under the same energy 
and centrality conditions constrained within a longitudinal rapidity window ly^ol < 0.5. In the 
right panels the abscissa for the proton and deuteron data are rescaled by a factor (indicated in 
the panel) allowing to superpose them on the triton data. At this beam energy this rescaling of 
the rapidity axes is the same for the rapidity distributions (which are in addition renormalized 
to the same area, lower right) and the flow (upper right). 

in Fig. |27]we show another instructive cut close to projectile rapidity, 0.85 < yo < 1.25. 
The sign of the mass hierarchy is inverted here. 

Attempts to extend 'simple' scalings of such data, as we showed in Fig. |25] for beam en- 
ergy 0.15y4 GeV, to higher beam energies, were moderately successful. In general, one can 
superimpose the V2{uto) functions of the three hydrogen isotopes if one uses momentum 
scaling factors in the ranges 1.5±0.1 (deuterons), resp. 2.0±0.2 (tritons). At close inspec- 
tion one finds that these scale factors vary somewhat within the indicated limits and, in 
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Fig. 26. Left: elliptic flow —^2(^*0) of hydrogen isotopes for Au+Au at various indicated incident 
energies with centrality 0.25 < 60 < 0.45 and in the rapidity window |yo| < 0.4. The smooth 
curves are three-parameter polynomial fits including even (2,4,6) exponent terms. 
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Fig. 27. Left: elliptic flow —V2{uto) of hydrogen isotopes for Au+Au at various indicated incident 
energies with centrality 0.25 < 6o < 0.45 and in the rapidity window 0.85 < |yo| < 1-25. 



contrast to the 0.15A GeV data, are larger than the factors needed to 'scale' the transverse 
rapidity spectra in a way similar to the case shown in the lower panels of Fig. [251 

Exploring the use of 'number scaling' (a nucleon number scaling here) in analogy to the 
(constituent) quark number scaling at RHIC [65], i.e. plotting V2/A versus the transverse 
kinetic energy per nucleon failed significantly and hence no such plots are presented here. 
As nucleons can be safely considered to be 'deconfined' at maximum pressure from the 
clusters they end up with at freeze out, one can as safely conclude that deconfinement 
does not necessarily lead to number scaling of observed elliptic flow. In the present energy 
regime from hydrodynamic scenarios (common flow velocity profile superimposed by ran- 
dom motion due to a common 'local' temperature) one does not expect number scaling, 
although a mass hierarchy is predicted [M]- As we shall see in the next section, microscopic 
transport theory is fairly successful in describing flow of various light clusters. 
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4-4 System size dependence of elliptic flow 



We study the system size dependence by comparing the rapidity dependence of the elhptic 
flow in ^°Ca + ^°Ca and Au+Au coUisions. This is shown in Fig.|28]for protons (left panels) 
and deuterons (right panels) at three different energies (top to bottom: OAA, l.OA and 
1.5A GeV) and centrality 0.25 < bo < 0.45. We observe a rather strong effect of the size of 
the system (our Ru+Ru data, not plotted here, are intermediate between the two systems). 
Not only the value at midrapidity is affected, but the whole shape of the curves: smooth 
flts of f 20 + f 22 Vo + ^24 iucluded in the figures. The heavier system is characterized by 

more 'compact' shapes. As we will see in section |5l this corresponds to the action of more 
repulsive mean fields. Using the stopping data from Ref. |40] as additional information, we 
can say that the larger densities reached in the heavier system, as implied by the larger 
stopping, generate elliptic fiow characteristic of a 'stiffer' system at compression time. 
As we shall see, the effect caused by the increased transparency of the lighter system is 
stronger than the difference between the so-called 'stiff' and 'soft' EoS options. 

4-5 Excitation functions for elliptic flow 

We have published earlier [29] an excitation function for the elliptic flow of Z = 1 particles. 
It has by now been joined up smoothly [16] to complementary lower energy data by the 
ALADIN-INDRA Collaboration. 

We show the earher FOPI data again here in the left panel of Fig. [29] (blue dots) together 
with our present evaluation which was done with a completely independent software and 
is seen to be in excellent agreement with [16]. The nominal impact parameter ranges are 
very similar: 6o=0.41-0.56 (for M3), 6o=0.45-0.55 (this work), both selected by multiplicity 
binning. Also, the same low momentum cut, Uto > 0.8 was applied in both evaluations. 
The smooth curve is a flt of a third degree polynomial in terms of log{E/u) to the present 
data that reproduces the V2 data with an average accuracy of 0.002. This comparison also 
conflrms that the quadrant method, section [2] and [39] is equivalent to taking < cos{2(j)) >, 
at least at SIS energies. However, the quadrant method could be used to design a way to 
correct for apparatus response asymmetries. As it turns out, due to the relative flatness 
of the f 2(1/0) dependences. Figs. [20], [23]and[28l the corrections to V2 at mid-rapidity, our 
present parameter ^20, are virtually negligible. 

With our present particle identifled data, and switching back to our sign change for positive 
'squeeze-outs', we can now plot (in the right panel of Fig. [29]) the separate excitation 
functions for masses one to four. To be able to join up the low and high energy data we 
again need to apply the common cut UtQ > 0.8. The smooth curves are polynomial flts 
along a logarithmic abscissa constrained to traverse zero at 4A GeV to be compatible with 
the proton data of ^72j . With this procedure, assuming for lack of better knowledge that 
the zero crossing at the high energy end is independent of particle type, we can summarize 
the complete range of energies where the phenomenon of 'squeeze-out' is present. Note that 
the various curves are not 'scalable' in a simple way along the — f2o ordinate: the locations 
of the maxima vary with the mass of the ejectile. This flnding could be influenced by the 
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Fig. 28. Elliptic flow —V2{yo) of protons (left) and deuterons (right) in ^'^Ca+^'^Ca (blue open 
triangles), and Au+Au (red dots) collisions for OAA (top), l.OA (middle) and 1.5A GeV (bottom) 
incident energy and centrality 0.25 < 6o < 0.45. A cut uto > 0.4 is applied. The smooth curves 
are fits of f 20 + V22 Ho + ^24 ut to the data. 
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Fig. 29. Excitation functions of elliptic flow in Au+Au collisions for centrality 0.25 < 60 < 0.45. 
Left panel: Elliptic flow V2 of Z=l nucleonic particles in Au+Au collisions as a function of incident 
energy. The present data (red open squares) are compared to the data published earlier |29] (blue 
dots). Right panel: Elliptic flow —V20 Protons: open black circles, deuterons: open red squares, 
mass 3: full green triangles, a's: pink asterisks. The smooth curves are polynomial fits along 
a logarithmic abscissa constrained to traverse zero at GeV [72]. A low pt cut was applied 
{uto > 0.8). This cut allows to join up the data from two different (low and high beam energy) 
runs. 

strong cut in uto which affects the heavier particles more strongly. More data between 2A 
GeV and 4A GeV would be useful: as shown before in this section and in [15] the onset of 
squeeze-out at the low energy end was found to be complex. 



4-6 Elliptic flow and isospin 



To study the isospin dependence of elliptic flow we proceed in analogy to our data analysis 
for directed flow and study either differences between the flow of '^H and ^He in a fixed 
N > Z system like Au+Au, Figs. |32] and |33l or look separately at various fixed ejectiles 
varying the isospin asymmetry of the system. Fig. [311 

However, before we proceed and, since the isospin effects are rather small, we show again 
in Fig. [30] the subtle influence of the corrections to the apparatus distortions, here for 
the V2{yo) functions. The left panel compares the flows of the isospin pair ^H/'^He after 
forward/backward symmetrization (see earlier. Fig. [1]), but ignoring the shift problems 
that were illustrated in Fig. [151 In the right panel the final data after full correction are 
shown. As we will see later in section [5l the small resulting differences over the entire 
rapidity interval are typical for a small increase of repulsive forces when switching from 
^He to '^H, a conclusion in line with our earlier observations using directed flow, section [31 

Similar to our observations for directed flow, we see no significant differences in the rapidity 
dependences f 2(2/0) for the various particles (protons, deuterons, tritons and ^He in Fig. 
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Fig. 30. Comparison of the elliptic flow —f 2(1/0) of (blue triangles) and '^He (red squares) in 
Au+Au collisions for l.OA GeV incident energy and centrality 0.25 < 60 < 0.45. A cut uto > 0.4 
is applied. Left panel: the data after applying the forward-backward symmetry constraint. Right 
panel: same data after additional ui-shift corrections. 

[3T|) when switching from ^^Zr + ^^Zr to ^^Ru + ^®Ru. The shown data are for 1.5A GeV 
incident energy with centrahty 0.25 < 60 < 0.45 and the general constraint Uto > 0.4. 

We see some interesting effects when we compare the isospin pair of mass three in Au + 
Au colhsions, Fig. [321 While the rapidity dependences for 0.4^4 GeV beams are still rather 
similar (upper left panel), a notable difference emerges at 1.5^4 GeV (upper right panel). 
The shape differences between and ^He data (for 0.25 < 60 < 0.45), larger \v2\ values at 
midrapidity for coupled to a more compact shape as function of rapidity, corresponds 
to effectively more repulsive field gradients. 

In order to quantize the differences into one parameter, dubbed V2n and shown in the lower 
two panels (left, an excitation function for fixed centrality, right, centrality dependence 
for fixed beam enegy), we have fitted the ^2(1/0) data in a limited range, \yo\ < 0.8, setting 
— "^2(2/0) = + V22 ■ Ho and defining V2n = V20 + |v22|- The parameter V2n increases both 
with the value and the curvature near mid-rapidity. With this parameter the difference 
between the two isotopes grows steadily with incident energy, being virtually zero at OA A 
GeV (lower left panel), but has the puzzling behaviour to disappear also at high energy 
(1.5^4 GeV) if the centrality is lowered (lower right panel). An interesting speculation 
at this point, explaining perhaps the large isospin difference at higher energies and larger 
impact parameters, could be that the effect is dominated by momentum rather than density 
dependence [73]. In [73] the simulations for Au+Au were hmited to 0.4^4 and 0.6^4 GeV and 
to our knowledge no covariant formalism was used, i.e. momentum dependences implied 
by Lorentz invariance were not included. 
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Fig. 31. Elliptic flow — 1^2(^0) of protons (upper left), deuterons (upper right), tritons (lower left) 
and ^He (lower right) in ^^Zr+^^Zr (blue open squares) and ^^Ru + ^^Ru (red dots) collisions 
for 1.5A GeV incident energy and centrality 0.25 < 60 < 0.45. A cut uto > 0.4 is applied. 



Using smooth fits to f 2(^*0) data similar to those presented in Fig. |26l we are able to 
show for the isospin pair in Fig. |33]the gradual evolution with energy starting from 0.12A 
GeV all the way to 1.5A GeV. Although the systematic error (0.1) of the original data 
for and ^He is rather large (and caused mainly by the uncertainty of the correction of 
symmetry violations by the apparatus), the evolution of the curves fitted to the data is 
surprisingly regular through-out this large span of energy and clearly shows the need for a 
systematic data taking. It will be an interesting task for the future to try to reproduce such 
data with transport codes dedicated to isospin dependences and capable of reproducing 
cluster yields. In view of the subtleties of such evolutions in isospin effects a confirmation 
of the present data would also be useful. 
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Fig. 32. Upper panels: elliptic flow —f 2(2/0) of (blue) and ^He (red) in Au+Au collisions 
for 0.4^4 and 1.5A GeV incident energies and centrality 0.25 < 60 < 0.45. Lower panels: elliptic 
flow characterizer V2n (see text) for (blue triangles) and ^He (red dots). Lower left: beam 
energy dependence for centrality 0.25 < &o < 0.45, lower right: centrality dependence for 1.5A 
GeV incident energy. 



4-7 Comparison to other data 



In Fig. [3l]we compare our data for Z = 1 and centrality 0.45 < &o < 0.55 to those of Ref. 
[50] which were the first 'squeeze-out' data to appear in the refereed literature. Originally 
these data were parameterized in terms of the 'number squeeze-out' parameter Rn which 
can easily be converted to the nowadays more common V2 using the simple relation 2v2 = 
(1 — Rn) /(! + Rn)- The authors of Ref. [SD] used the flow reference frame mentioned in the 
beginning of this section, but to our knowledge no reaction plane correction was applied. 
Therefore, to compare, we have omitted the correction as well ( for the corrected data, see 
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Fig. 33. Elliptic flow —f 2(1^40) of "^H (blue) and ^He (red) in Au+Au collisions for different 
indicated incident energies and centrality 0.25 < 60 < 0.45. The smooth curves were obtained 
from three parameter fits of even polynomials (excluding the constant term) to the data. The 
low momentum parts (below uto = 0.6 for the upper three panels, below 0.4 for the other panels) 
are extrapolations of the data. 

our Fig. |29|) . Rotation into the flow axis should maximize the elliptic flow (see also |15]) 
and indeed, for incident energies below 0.5 A GeV our non-rotated data show signiflcantly 
less elliptic flow values. However, at higher energies our data (despite not including the 
stronger flow of Z = 2 fragments) suggest increasingly more flow than the rotated data. 
This is not understood. 

In Fig. [25] we compare to proton flow data from the KaoS Collaboration [71]. These 
data, besides the kaon data [33'J34J, played an important role in some past theoretical 
efforts to constrain the EoS j23p27j . We expect the KaoS data to have a better momentum 
resolution, but the equally important centrality selection and reaction plane resolution 
should be superior in our large acceptance apparatus. 

At an incident energy of 0.4A GeV we have unfortunately no direct match of centralities 



43 




0.4 0.8 1.2 
beam energy (A GeV) 



1.6 



Fig. 34. Elliptic flow data for Au+Au collisions 
as function of beam energy. Open red squares: 
data from Ref. [50] for Z = \ and 2 particles 
in the flow reference frame. Closed black dots: 
present data in the beam axis reference frame 
for Z = 1 and centrality 0.45 < 6o < 0.55. Both 
types of data are uncorrected for reaction plane 
resolution. 




Fig. 35. Elliptic flow —^2(^*0) of protons in Au+Au collisions at 0.4^ GeV (left panel) and l.OA 
GeV. The present data {MUL selected) are compared to the data adapted from Ref. [74j (KaoS 
coll.). The two panels show results for different indicated average centralities. The KaoS data 
correspond to blue filled squares. The smooth curves to the FOPI data are fits with a polynomial 
Vt22 • Uto + 'Vt23 • Uto + '^i24 " u^q. In the right panel the FOPI fit is seen to also describe the KaoS 
data if it is downshifted by 0.03. 



and show therefore the KaoS data for bo ~ 0.45 together with our data for two (average) 
centralities 60 = 0.36 and Bq = 0.50. For large enough momenta there is compatibility of 
the two sets of data, but we do not find evidence for a sign change at low uto- The 1.0 A 
data (right panel) match almost perfectly in centrality, but (see the caption) the KaoS 
data are systematically downshifted from the FOPI data and again suggest a sign change 
at low momenta not seen by us. 
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5 Transport model simulations 



Before showing simulations of our data with a transport code, we wish to stress in Fig. [36] 
the importance of having isotope separated data to do so. The figure shows for Au on Au 
collisions the various direct and elliptic flow projections that we have studied in sections |3] 
and m for protons (black open circles) and for the sum of the three hydrogen isotopes. The 
differences are quite significant and complex. Clearly, if the relative yield of the various 
isotopes is not predicted correctly in the complete three-dimensional momentum space the 
interpretation of data that are only charge separated is not reliable. It is also important to 
realize this when comparing proton or neutron flow and hydrogen flow to learn something 
about isospin dependences. 

As transport code we use a quantum molecular dynamics code, IQMD. The IQMD ver- 
sion we use in the present work largely corresponds to the description given in Ref. |12] . 
More specifically, the following 'standard' parameters were used throughout: L = 8.66 fm^ 
(wave packet width parameter), t = 200 fm (total propagation time), K = 200,380 MeV 
(compressibility of the momentum dependent soft, resp. stiff EoS), Egym = 25p/po MeV 
(symmetry energy, with po the saturation value of the nuclear density p). The versions 
with K = 200, resp. 380 MeV are called IQMD-SM, resp. IQMD-HM. The clusterization 
was determined from a separate routine using the minimum spanning tree method in con- 
figuration space with a clusterization radius Rc = 3 fm. 

We are aware that in the wake of intense nuclear symmetry energy research more recently 
'improved' transport codes have been developed that address isospin dependences by in- 
creasing the relevant parametric flexibility but not necessarily the general consensus what 
to conclude: ImQMD [80], IBUU04 [21], UrQMD [15]. Our motivation to stick for now 
to the present IQMD version is that we have used it in two of our earlier papers on pion 
systematics [3Dj and the systematics of central collisions [40j and want here to test IQMD 
knowing how it performed on the observables published earlier. We believe that convinc- 
ing conclusions on basic nuclear properties require a successful simulation of the full set 
of experimental observables with the same code using the same physical and technical 
parameters. 

The two options of purely phenomenological cold nuclear EoS that we use are plotted 
in Fig. [37] and confronted with a 'microscopic' (Dirac-Brueckner-Hartree-Fock, DBHF) 
calculation [75] for symmetric matter. Also included from the same theoretical work is 
the cold EoS for pure neutron matter. These theoretical predictions are fully in-line with 
what is known so far on both EoS in the vicinity of the saturation density po. It is seen 
that in the density range relevant for SIS energies (up to p/po = 2.5) our 'soft' version, 
SM, is rather close to the theoretical calculation. In this density range SM is also close 
to the EoS dubbed AP4 in Refs. [76][77] which is due to a calculation [78] based on the 
Argonne vis two-body potential combined with the Urbana IX three-body potential and 
relativistic boost corrections (it is dubbed Al8+(5f-|-UlX in Ref. [78]). In a recent astro- 
physical measurement [77j it was shown that AP4 was compatible with the discovery of 
a two-solar-mass neutron star. In order to support such a massive object the nuclear EoS 
has to be sufficiently stiff. (Actually AP4 has to be shifted down by about three MeV to 
be in good accord with the binding energy per nucleon at saturation density po). 
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Fig. 36. Comparison of hydrogen and proton flow in 0.4^4 GeV Au+Au coflisions with centrahty 
0.25 < 6o < 0.45. Upper left panel: rapidity dependence of vi with a constraint uto > 0.4. Upper 
right panel: uto dependence of vi with a constraint 0.4 < |yo| < 0.8. Lower left panel: rapidity 
dependence of —V2 with a constraint uto > 0.4. Lower right panel: uto dependence of —V2 with a 
constraint |yo| < 0.4. 

The EoS for intermediate values (0 to 1) of the asymmetry parameter 6 = {pn — Pp) /p is 
expected in most of the literature to have a nearly quadratic dependence on 6, E{p, 6) = 
Eo{p) + Easyip) ■ 5^, a relation confirmed for example in [75]. The 5^ values available in 
earthly laboratories are rather modest: 5^ = 0.0447 for the readily available 208p]^_j_208p]^ 
and 5^ = 0.0478 for a popular radioactive beam system in a foreseeable future, ^"^^Sn+^^^Sn. 
The binding energy curves for such systems are much closer to symmetric nuclear matter 
than to neutron matter. 

We have sorted our simulations by incident energy, as the physics, and the associated 
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Fig. 37. The two equations of state HM (red dashed) and SM (blue continuous) are compared 
to symmetric (black crosses) and neutron matter (blue circles) EoS obtained using DBHF theory 
(after [75j). 

problems understanding it, vary with this externally diallable parameter. 
5.1 Comparisons at 0.15 A GeV. 

We start at an incident energy still close to the balance energy, E = 0.15 A GeV. In Fig. 
l3S]we compare the rapidity dependences of elliptic flow in Au + Au collisions of centrality 
0.25 < 60 < 0.45. Proton data are shown in the left panel, deuteron flow is plotted in 
the right panel. As for all our low energy data, a constraint Uto > 0.8 is applied which 
is expected to influence the results significantly, but, as mentioned earlier, preserves the 
symmetry around mid-rapidity. The constraint is of course also applied to the simulation. 
In this and all following simulations we always use ERAT to select centrality both in 
the experimental and theoretical events. In the figure both the HM and the SM options 
are shown together with the data (black dots with systematic errors). To improve on 
statistical significance we have fitted two-parameter parabola's f 20 + "^22 ■ Vo to the various 
distributions. At this incident energy such shapes spanning beyond the full scaled rapidity 
gap (-1 to +1) are seen both in the experiment and the simulation. 

Some remarks can be made: 
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Fig. 38. Rapidity dependence of elliptic flow —V2 of protons (left) and deuterons (right) in 
Au+Au collisions at E = 0.15 A GeV for 0.25 < 6o < 0.45. Transverse 4- velocities uto below 0.8 
are cut off. The curves are fits of V20 + f 22 2/0 *o data. The black dots are the experimental 
data. The error bars of the simulated data are statistical assuming the fitted 2-parameter shape 
is correct and for clarity are indicated only for the IQMD-SM calculation (blue full curve). The 
red dashed curve represents the IQMD-HM calculation. 

1) At this low energy and with the relatively high cut on Uto the difference between the 
various EoS predictions is rather modest. 

2) For protons the correct emergence of a positive squeeze-out, i.e. a dominant repulsion, 
is missed and probably influenced by the effective range of nucleon forces used, that in 
IQMD depends on partly technical parameters such as the 'wave packet width' parameter 

3) Deuteron, and more generally cluster flow gives more prominent signals and, surpris- 
ingly, is better reproduced than proton flow. 



5.2 Comparisons at 0.25 A GeV. 



Moving to 0.25y4 GeV where definitely repulsive behaviour is established in the Au + 
Au system, we show in Fig. [39] directed flow data for protons, deuterons, tritons and 
^He together with a simulation using the SM option. For the chosen centrality, 0.25 < 
bo < 0.45 and cutting off momenta corresponding to Uto < 0.8, we find an amazing 
agreement between experiment and simulation (note the changing ordinate scales in the 
various panels). In view of the fact that absolute yields of such clusters are not well 
reproduced by the present IQMD version [18], this is not trivial. Of course, as mentioned 
earlier, missing the relative yields, the flows of Z = 1 or 2 fragments, not separated by 
mass, would be off. 

The success of IQMD at incident energies below 0.4^4 GeV where the sensitivity to the 
EoS is still weak can be seen as an important step towards making sure that some of the 
technical parameters of the code have appropriate values. The possible influence of these 
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Fig. 39. Directed flow vi of protons, deuterons (upper panels) and of and "^He clusters emitted 
in Au+Au collisions at E/u = 0.25^ GeV and centrality 0.25 < 6o < 0.45. Transverse 4- velocities 
uto below 0.8 are cut off. The small filled black circular dots represent the experimental data. 
The smooth curves are least squares fits of vi = vu ■ yo + V13 ■ Uq to the simulation results (open 
blue triangles) with a soft EoS. 



parameters was discussed in [12] and more recently, for low E/u, in [T5] . 



5.3 Comparisons 0^0.4^4 GeV 



At an energy of OAA GeV we start considering more seriously the differences between 
soft and stiff EoS options and use for that the extended information of both directed and 
elliptic flow, watching the rapidity as well as the transverse momentum dependences. 

In the three panels of Fig. HOlthis is done for protons, in Fig. |4T]it is shown for deuterons 
and in Fig. |12] for tritons. Again the centrality is 0.25 < 60 < 0.45, but at this higher 
energy we can use the less restrictive constraint Uto > 0.4 for the momentum integrated 
rapidity dependences a fact which tends to increase the sensitivity to the 'cold' EoS: low 
momentum particles are more affected by mean fields. The momentum dependences (upper 
right panels) are for the rapidity interval |?/o| < 0.4. In all the panels the smooth curves 
are polynomial least squares fits using the same fitting parameters (at most three) and 
abscissa ranges for data and simulation. As shown specifically for the f 2(^*0) data such 
fits reproduce the experimental data well within their error bars. 
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Fig. 40. Flow of protons in Au+Au collisions at E = OAA GeV and centrality 0.25 < bo < 0.45. 
Smooth fits to simulation data (red dashed: stiff EoS HM, blue full: soft EoS SM) are compared to 
FOPI data (black filled dots). Transverse 4- velocities uto below 0.4 are cut off. The uto dependence 
of {—V2) (upper right panel) is shown within a cut |yo| < 0.4. 

Looking at the figures we can note the following: 

1) There is a clear preference for the soft version, SM, of the EoS. 

2) The conclusion is the same for elliptic flow as for directed flow. 

3) Amazingly, again, the rendering of the proton data is less well achieved than for the 
deuteron and it is best for triton clusters. 

4) Looking especially at V2{yo) we see an increasing sensitivity to the EoS as the mass of 
the ejectile is increased (right lower panels). 

While the comparison of experimental data with the simulation, in particular the SM 
version of it, so far looks quite encouraging, the attentive reader will notice that we have 
not, in the last three figures, compared to the momentum dependence of directed fiow. 
This is done for the three hydrogen isotopes in Fig. |13] for the SM option showing an 
intriguing discrepancy with experiment that is growing as the emitted particle becomes 
lighter. For a centrality 0.25 < 60 < 0.45 we show the projection of vi on the Uto axis in a 
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Fig. 41. Flow of deuterons in Au+Au collisions at E = OAA GeV and centrality 0.25 < bQ < 0.45. 
Smooth fits to simulation data (red: stiff' EoS HM, blue: soft EoS SM) are compared to FOPI 
data (black filled dots). Transverse 4- velocities uto below 0.4 are cut ofi^. The uto dependence of 
(— ^2) (upper right panel) is shown within a cut |yo| < 0.4. 



rapidity interval 0.3 < \yo\ < 0.7. While for the data the slope of vi{uto) at low momenta 
varies significantly with the mass of the particle, being steepest for the triton, this effect 
is completely missing in the simulation. In particular the gentle slope of Vi for protons is 
not reproduced. 

Searching for reasons why the proton flow data are more poorly reproduced by the sim- 
ulations, it is instructive to come back to Fig. [36] which compared measured hydrogen 
{Z = 1) flow with the measured proton {A = 1) flow: one sees that the differences between 
[Z = 1) and [A = 1) look qualitatively very much like the differences in Figs. 140) and 143) 
between the SM calculation and the proton data. This suggests that the lack of sufficient 
clusterization (rather than an inappropriate EoS) in the simulation is at the origin of the 
deficiency of the theory. The deuteron and the triton data are less influenced by this effect 
and hence better reproduced. 
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Fig. 42. Flow of tritons in Au+Au collisions at E = OAA GeV and centrality 0.25 < 60 < 0.45. 
Smooth fits to simulation data (red dashed: stiff EoS HM, blue full: soft EoS SM) are compared to 
FOPI data (black filled dots). Transverse 4- velocities uto below 0.4 are cut off. The uto dependence 
of {—V2) (upper right panel) is shown within a cut |yo| < 0.4. 



In view of Fig. Il3] another interpretation cannot be excluded as well: It could be that 
an important quantum mechanical effect is missing at freeze-out in the (quasi-classical) 
simulation: the Fermi-motion of the clusters in the medium, which is expected to be 
largest for single nucleons. As Fermi motion from a flow (i.e. hydrodynamic) point of 
view is chaotic, it will tend to 'soften' the visible flow just as thermal motion does. As a 
consequence f 1(1/0) for protons (Fig. HQ]) is more affected than fi(?/o) for tritons (Fig. H2l): 
Brownian motion is less imporant for heavier clusters. These conjectures require further 
investigation as they might shed additional light on the clusterization mechanism. As we 
shall see, this deficiency of the simulation is not unique to 0.4^4 GeV. Also, if one looks 
back at Figs. [39] to Il2] one notices that the f i(?/o) data of protons and deuterons are better 
described by the simulation for 0.25A GeV than for OAA GeV. With the information from 
Fig. Us] we can say that this is probably due to the higher cut on uto operated on the 0.25^4 
GeV data (0.8 instead of 0.4). 
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Fig. 43. Flow vi{uto) of hydrogen isotopes in Au+Au collisions at E = OAA GeV and centrality 
0.25 < 60 < 0.45. The rapidity is constrained to the interval 0.3 < |yo| < 0.7. The experimental 
data (full symbols) are compared to simulations using IQMD-SM. Smooth fits to FOPI data are 
used to extrapolate to uto = 0. 

5.4 Comparisons at 1.0 A GeV 



We proceed at an energy of lA GeV in a similar way to what we did at 0.4^1 GeV showing 
flow data for protons, Fig. HH and deuterons, Fig.|45]and also the remarkable disagreement 
with the simulation concerning the mass dependence of vi{uto), Fig. |46l The same cuts in 
scaled units as we used at OAA GeV are applied. 

As before, the SM option gives a better description of the data, especially of the momentum 
integrated rapidity dependences. For the deuteron data, due to the higher sensitivity of 
clusters, the preference for SM over HM is even more convincing. In the upper left panels 
of Figs. SH and SS] we have also plotted the rapidity distributions of the yields, applying 
here also the Uto > OA constraint both for the fi(|/o) and ^2(2/0) data. Comparing to the 
simulation, we see again the now well known overestimation of single nucleon emissions 
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Fig. 44. Flow of protons in Au+Au collisions at E = l.OA GeV and centrality 0.25 < 6o < 0.45. 
The simulation data from IQMD-SM are plotted as smooth full curves (blue), from IQMD-HM as 
dashed curves (red) FOPI data as black dots. The upper right panel shows the uto dependence in 
the indicated yo bin, the lower two the rapidity dependences integrated over uto, but constrained 
to Uto > 0.4. Directed flow, vi, is on the left, elliptic flow, —V2, on the right. The rapidity 
distributions in the upper left panel are also constrained to uto > 0.4. 



accompanied for deuterons by the underestimation of deuteron yields. The more interesting 
point, made by us earlier [^HPO] , is that when the system is stiffer (HM, dashed curves) 
this effect is even more pronounced: there are less clusters. We have associated this with 
the fact that stiff systems achieving less maximal compression during the collision, undergo 
later a less efficient cooling in the expansion period. 

In contrast to the rapidity dependences, the momentum dependences of proton flow are 
not so well reproduced. In Fig. |44]for example it is seen that for high momenta {uto > 1) 
the preference for a given EoS is no longer clear. In Fig. |16], where we show vi{uto) for the 
three hydrogen isotopes in the rapidity interval 0.3 < \yo\ < 0.7, we are stressing again 
the qualitative discrepancy between experiment (left panel) and simulation (right panel) 
by rescaling the Vi of deuterons and tritons so that they agree with the proton value at 
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Fig. 45. Flow of deuterons in Au+Au collisions at E = l.OA GeV and centrality 0.25 < 6o < 0.45. 
The simulation data from IQMD-SM are plotted as smooth full curves (blue), from IQMD-HM as 
dashed curves (red) FOPI data as black dots. The upper right panel shows the uto dependence in 
the indicated yo bin, the lower two the rapidity dependences integrated over uto, but constrained 
to Uto > 0.4. Directed flow, vi, is on the left, elliptic flow, —V2, on the right. The rapidity 
distributions in the upper left panel are also constrained to uto > 0.4. 

high scaled momenta. (The factors are indicated in the figures). While the IQMD-SM 
predictons are then superimposable for the three isotopes, this is clearly not the case for 
the experimental data. The latter show gradually steeper slopes at low uto when shifting 
from protons to deuterons and then to tritons. 

5.5 Comparisons at 1.5 A GeV 

We come finally to 1.5 A GeV, the highest incident energy currently available at SIS for the 
heavy system Au+Au. For the same centrality as before, we show the four flow projections 
that we have studied before. With the exception again of the Vi{uto) functions (around 
Ho = 0.5) and the high momentum part of f 2(^*0) for protons, the description of the proton 
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Fig. 46. Flow vi{uto) of hydrogen isotopes in Au+Au collisions at E = l.OA GeV and centrality 
0.25 < 6o < 0.45. The rapidity is constrained to the interval 0.3 < |yo| < 0.7. The experimental 
data (left panel) are compared to simulations using IQMD-SM (right panel). The ordinates of the 
deuteron and triton data are rescaled as indicated in the panels. Note the rescaling is different 
in the simulation. Smooth fits to FOPI data are used to extrapolate to uto = 0. 

data, Fig. |47l and the deuteron data, Fig. HHl is close to perfect if we use the SM version. 
Looking, specifically, at the elliptic flow f 2(2/0) for deuterons, we note that IQMD-SM is 
doing an excellent job from 0.15 A GeV (Fig. [3H1) to 1.5 A GeV (Fig. HHl) Au on Au, i.e. 
over an order-of-magnitude change in energy, a remarkable achievement of the transport 
code used. We can say this despite the 'failures' concerning the yields of clusters and the 
mass hierarchy of the vi{uto). So more work is needed, but an encouraging start is done. 

5. 6 Excitation functions 

In Figs. [38] to HHl we have tried to reproduce with simulations the full spectral information 
contained in the various projections of directed and elliptic flow. This is the most convinc- 
ing way of checking the adequacy of the simulation. To be able to represent such data in a 
readable form we had for a given plot to fix many parameters, such as system size, incident 
energy, centrality, particle type and finally the projection constraints. To get a somewhat 
better overview of more general trends, such as excitation functions, it is desirable to try 
to characterize flow with just one parameter. For elliptic flow this is generally done by 
just showing its value at mid-rapidity as we have done in Fig. [291 For directed flow the 
mid-rapidity slope dvi{yO) / di/Q = vn or its momentum- weighted alternative [281179] . u^oi 
in our scaled notation, is introduced. Excitation functions for vn were shown in Fig. [T4l 
While this more compact form of displaying flow data is useful to display general trends, 
one has to be aware that 

a) information is lost (we showed the instructive shape changes of f 2(^/0) associated with 
changes in the stiffness of the assumed EoS), and 
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Fig. 47. Proton flow for centrality 0.25 < 6o < 0.45 in collisions of Au+Au at 1.5A GeV. The 
left panels show directed flow, elliptic flow is plotted in the right panels. The lower panels show 
the rapidty dependence, while the upper panels illustrate the momentum dependence. Both the 
simulations (SM and HM) are shown together with the data (black dots). For further details see 
text. 



b) a less 'robust' comparison is made: the midrapidity directed flow slopes depend on some 
technicalities such as the range of the fltted data and the degree of the polynomials fltted, 
and clearly 'local' systematic errors become more important when narrower data samples 
are taken. 



With this proviso in mind we present in Fig. |19] a systematics of slopes Uxoi for deuterons 
as function of incident energy varying the centrality as indicated in the various panels. 
We have used the same procedure throughout: fltting the interval —1.1 < < 1.1 with a 
polynomial including a linear and a cubic term. We choose deuterons rather than protons 
for two reasons: 

1) as shown before, the flow signal is larger, and 

2) in this energy range deuterons are less likely than protons to be resulting in the last 
stages of the expansion from either late decays of larger clusters or decay of baryonic res- 
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Fig. 48. Deuteron flow for centrality 0.25 < 6o < 0.45 in collisions of Au+Au at 1.5A GeV. The 
left panels show directed flow, elliptic flow is plotted in the right panels. The lower panels show 
the rapidty dependence, while the upper panels illustrate the momentum dependence. Both the 
simulations (SM and HM) are shown together with the data (black dots). For further details see 
text. 



onances. 

As shown earlier, the energy range we choose (0.4 — X.'bA GeV) for this figure is character- 
ized by rather small changes of flow if scaled units are used. We ignore in the figure therefore 
small variations with incident energy and draw straight and constant lines through the 
data (and the HM and SM calculations) that characterize the average behaviour. Since in 
the previous figures we have shown comparisons almost exclusively for one fixed centrality, 
0.25 < 6o < 0.45, we focus our interest here on varying the centrality. The upper right 
panel confirms the almost perfect agreement of the data with the SM option. While this 
agreement can be extended to a more peripheral interval, 0.45 < 6o < 0.55, in the upper 
left panel, the data shown in the lower panels seem to suggest an even 'softer' behaviour 
than SM for the most central collisions. This strengthens the arguments tending to discard 
a stiff option like HM, but also confirms that more work is necessary to still better under- 
stand the 'details'. In particular it must be checked if stopping is described properly (see 
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Fig. 49. Excitation function for midrapidity slopes Uxoi of directed flow of deuterons in collisions 
of Au+Au with centrality 0.45 < 5o < 0.55 (upper left panel), 0.25 < 6o < 0.45 (upper right 
panel), bo < 0.25 (lower left panel) and bo < 0.15 (lower right panel). A low-momentum cut-off 
was done {uto > 0.4). The FOPI data (black full circles) are compared to simulations with 
IQMD-HM (red open squares) and IQMD-SM (blue open triangles). The straight lines represent 
averages over the shown energy range. 

later). Also the missing Fermi motion in the simulation, discussed as possible reason for 
the incorrect rendering of vi{uto), especially of protons. Figs. Il3]and|l6l probably causes 
part of the overestimation of the Vu slopes even when the SM option is used. As noted 
before, interpreted classically, Fermi motion or Brownian motion, has the effect of higher 
apparent temperature lowering the slopes Vu in the experiment. 

We have pointed out before the relatively small variation of directed flow with the cen- 
trality bo and associated this tentatively with the higher stopping in very central collisions 
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Fig. 50. Midrapidity slopes of directed flow of hydrogen isotopes in collisions of Au+Au 
with centrality < 0.15 at six different indicated beam energies. A low-momentum cut-off was 
done {uto > 0.4). The FOPI data (black dots) are compared to simulations with IQMD-HM 
(red open squares) and IQMD-SM (blue open triangles). The straight lines represent linear fits 
guiding the eye. 
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compensating the loss of geometrical asymmetry. Also we saw that heavier clusters showed 
more flow and this flow was more sensitive to the EoS than flow of single nucleons. Us- 
ing the more compact form of flow description, we confirm these trends in Fig. |50] which 
presents a midrapidity slope systematics (uxoi) for the most central collisions (60 < 0.15, 
2% sharp-radius cross sections) as a function of ejectile mass for six indicated incident 
energies. For such pronounced centrality, using the ERAT observable, it is particularly 
important that the simulations are done using the same selection methods as in the exper- 
iment. Looking first at the lowest right panel, which shows Au on Au data for the highest 
energy, 1.5A GeV, we note that HM predictions for tritons are almost 100% (!) too high: 
this is a large signal against HM. The preference for SM is true for all E/u, but at close 
inspection one sees some fluctuations comparing the straight lines drawn through the SM, 
HM and experimental data. These fluctuations seem to be of numerical origin and call for 
further efforts. 

5. 7 The other observables: radial flow and stopping 

To raise the confidence in the physical adequacy of the simulation, it is important not to 
limit the confrontation with experiment to directed and elliptic flow. Besides the chemical 
composition at freeze-out [IQ] there is another related observable, the radial flow which 
is the main signal of the expansion-cooling process. Even more important is the so-called 
stopping. A rather detailed presentation of these observables has been done earlier [ID]. 

In Fig. |5l] we present a systematics of radial flow for six incident energies together with 
transport code simulations. What is plotted is the scaled variance, varxrriQ of constrained 
(to \yo\ < 0.5) transverse rapidity distibutions |10] as function of mass for the hydrogen 
isotopes. Note that we are looking at the most central collisions here as for the previous 
(and the next) figure. In a globally thermal model (no flow, complete equilibration) these 
variances should be (approximately, ignoring relativistic effects) proportional to A~^, the 
inverse of the ejectile masses. So then, when multiplied with the mass A, as done in the 
figure, they should not depend on mass. The fact that they grow with mass is interpreted 
as 'radial flow' in most of the literature. Such radial flow is also predicted by the IQMD 
calculation as can be seen in the figure and, obviously, depends very little on the stiff- 
ness of the EoS. This was known before, but is seen here to reproduce, without any ad 
hoc adjustment, the data very well at the highest energies (the two lower panels) and 
overestimating the mass trend somewhat at the lowest energies (upper two panels). For 
hydrodynamically oriented people we add a word of caution: due to the observed stopping 
hierarchy as function of mass (see Fig. [52]) the interpretation of slopes as seen in Fig. [51] 
exclusively in terms of a common 'radial-flow' velocity (profile) cannot be quantitatively 
correct. 

One can safely expect that maximal density is only achieved if full stopping has occurred. 
On the other hand incomplete stopping can be misinterpreted as 'ultra-soft' EoS since 
it is associated with less static pressure. Stopping must therefore be studied with great 
care. Our registered events, when sorted according to ERAT are in fact sorted according 
to the degree of stopping. While ERAT is a global event qualifier affected partially by 
apparatus effects (like multiplicity it is only used for selection), we defined in our earlier 
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Fig. 51. Mass-weighted scaled variances varxiriQ of constrained transverse mid-rapidity distribu- 
tions for hydrogen isotopes in coUisions of Au+Au with centrahty 6o < 0.15 at six different indi- 
cated beam energies. The FOPI data (black dots) are compared to simulations with IQMD-HM 
(red open squares) and IQMD-SM (blue open triangles). The straight lines represent linear fits 
to the simulated data. 
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Fig. 52. Stopping varxz for hydrogen isotopes in collisions of Au+Au with centrality bo < 0.15 at 
six different indicated beam energies. The FOPI data (black dots with error bars) are compared 
to simulations with IQMD-HM (red open squares) and IQMD-SM (blue open triangles). The 
straight lines represent linear fits to the simulated data guiding the eye. 
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work [lO] an observable varxz=varx/varz where varx is the variance of the transverse 
1-dimensional rapidity distribution (in an arbitrary direction x) and varz holds for the 
rapidity distribution in the beam direction (the commonly used longitudinal distribution) . 
varxz is determined individually for each particle type and corrected for An acceptance. It 
was observed PU] to follow a 'stopping hierarchy': heavier clusters were associated with 
higher partial transparency (smaller varxz). Our simulations also predict this hierarchical 
trend for the hydrogen isotopes as evident from Fig. [521 In contrast to radial flow, there is 
some sensitivity to the EoS [IQ] as documented here also in Fig. |52l It is intuitively clear 
that a stiff repulsive mean fleld will favour stopping. We note in this context that our 
version of IQMD does not consider in-medium modiflcations (besides Pauli blocking) of 
short-range nucleon-nucleon interactions and cross sections. The strong axially symmetric 
squeeze-out {varxz > 1) predicted very early [3j from ideal hydrodynamics is not observed: 
it is hidden by (forward peaked !) surface transparency effects that were not considered at 
the time: as one can see from the flgure the varxz values predicted by microscopic transport 
are well below one in agreement with experiment (black dots in the flgure). On the other 
hand the hierarchy of stopping is somewhat less pronounced in the experimental data and 
in particular the calculation misses the absence of the hierarchy at 0.4 — 0.6^4 GeV where 
stopping was found |66] to be maximal. 



6 Conclusion 

The present work complements a systematic study of heavy ion collisions at SIS energies 
(0.09 — 1.5A GeV) with the aim of establishing a broad data base for theoretical efforts to 
deduce constraints on the nuclear EoS and other nuclear properties, such as viscosity. In 
the frame of this systematics we have published [39] extensive data on pion production, 
pion spectra and pion flow. More recently, we have presented for central collisions an 
overview of the chemical composition of the expanding flreballs, a systematics of stopping 
and the evidence for radial flow in [ID] . Here we have presented a systematics of azimuthally 
asymmetric fragment emissions in terms of the well established parameters vi, directed 
flow, and V2, elliptic flow. As already pointed out in some of our earlier work, [SUSS], both 
vi and V2 are complex multidimensional functions: as a consequence we found it necessary 
to extract from our data dependences on incident energy, system mass and charge, ejectile 
mass and charge, centrality and, of course, longitudinal rapidity and transverse momentum. 

While the energy range spans more than an order of magnitude change, the system sizes 
were varied from ^"Ca+^^Ca to ^^^Au-|-^^^Au and the isospin of the system was separately 
varied using ^^Ru+^^Ru and ^^Zr-|-^^Zr. The centrality selection method was based pri- 
marily on binning the distribution of the ratio of transverse and longitudinal energies, 
ERAT, i.e. on the degree of stopping achieved. We stress that we use well deflned cross 
section intervals, corresponding to scaled impact parameter ranges, bo < 0.15, bo < 0.25, 
0.25 < &o < 0.45, and 0.45 < &o < 0.55. We cover a large, sharply deflned two-dimensional, 
ivo, Uto), phase space both for vi and V2. AH these features should ease precise duplication 
by transport model simulations, a necessary condition to extract reliable nuclear matter 
properties from the experimental data. Of particular interest are the rapid rise of flow 
below O.SA GeV, followed by a relatively flat behavior above. It will also be instructive to 
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try to reproduce quantitatively the correlations between system size dependent flow and 
stopping [i0|66] . 



Using as transport code IQMD, the same code/version that we used in our earlier studies 
on different observables, we have done an orientational study to see how well we could 
reproduce the flow data and to which degree they were sensitive to the assumed EoS. To 
our surprise, although the used version of IQMD does not well reproduce cluster yield data, 
it reproduces well the flow of light, mass and charge identifled clusters, namely deuterons 
if a soft equation of state is assumed. As a matter of fact, our confrontation with data 
virtually excludes a stiff EoS. This strengthens the case of a soft EoS proposed in earlier 
works |31f35f36] using different observables. Another gratifying result is that we do not 
have any difficulty reproducing both directed and elliptic flow with the same assumptions. 
At the highest energy studied here, we were also able to reproduce well the radial flow 
and the stopping for the various hydrogen isotopes. We note that our version of IQMD 
did not assume any in-medium modiflcations of the microscopic cross sections other than 
Pauli blocking. No ad hoc modiflcation to achieve agreement with our data was done. 

More speciflcally, for elliptic flow, V2{yo), we were able to check the influence of the EoS on 
the shape of the function in a large rapidity range, not just at mid-rapidity. This allowed us 
to state that the differences of ^2(2/0) between and '^He are consistent over the complete 
rapidity gap with a more repulsive force acting on the neutron richer cluster. Our directed 
flow data on the mass three clusters are also consistent with this conclusion. 

Still, more work needs to be done. On the experimental side a conflrmation of the isospin 
effect will be helpful. The future international facility FAIR at GSI, Darmstadt should 
allow to access higher nuclear densities than has been possible here. Theoretically it would 
be desirable to reproduce the cluster yields without loosing the good performance on flow 
observables. Also, using internally consistent theories such as [T9|, a further flne tuning of 
our stopping data is needed. More subtle effects, to be further investigated, are possible 
correlations of flow to clusterization visible in the momentum dependences of single nucleon 
flow and, also, the experimental evidence of quantum effects on flow, such as Fermi motion 
subsisting at freeze-out time. Fermi motion tends to get 'lost' in a semiclassical transport 
calculation. Such deflciency leads to overestimation of the predicted flow. All of this is an 
ambitious task beyond the scope of the present experimental work. Hopefully there will be 
convergence of the various existing transport codes in the future. An exciting prospect is 
a possible joining up of information from heavy ion data with astrophysical observations 
on neutron stars [ST|52] . 
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Appendix 



The correction factors [52] given in the table are defined by 

vi = v[/di and V2 = v^Ci 

where the primed quantities are the uncorrected flow parameters and the index i distin- 
guishes the various i?i?yl T selected centrahties, in order 
bo < 0.15, bo < 0.25, 0.25 < 6o < 0.45 and 0.45 < 6o < 0.55. 
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